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KINETIC PERFORMANCE HANDBOOK 


"by 

L. S. Bender, V. J. Sarli, and W. G. Burwell 


ABSTRACT 


Presented in this handbook, prepared "by the United Aircraft Research Labora- 
tories under Contract NAS 3-11225, are self-consistent kinetic performance and 
aerodynamic degradation charts to permit rapid evaluation of the deliverable per- 
formance for five liquid propellant combinations, H2-F2, Aerozine 5O-N2OI1., 

CH^-Flox (82.6$ F2), and B2Hg-0F2 over wide ranges of operating conditions and 
for different nozzle geometries and engine sizes. The conditions generally en- 
compass a range of oxidizer/fuel ratios surrounding the anticipated optimum mixture 
ratio, and include combustion chamber pressures from 100 to 1000 psia (50 to 500 
psia for B2Hg-0F2)* Charts are presented for both conical and bell nozzle contours 
for values of nozzle exit area ratio between 20 and 100 and for nominal thrust 
levels between 100 and 1,000,000 pounds. Also included is sufficient information 
for the user of this handbook to estimate the effect on kinetic performance loss 
of variations in reaction rate constants relative to those used in the handbook. 

The procedures and reaction rate constants employed herein are consistent with 
those recommended by the ICRPG Performance Standarization Working Group. 
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INTRODUCTION 


Extensive effort has been expended in recent years to assess "real gas" flow 
properties in subsonic-supersonic expansion nozzles. A particularly large amount 
of work has been performed with the objective of establishing chemical nonequili- 
brium flow effects in such nozzles, a prerequisite to predicting, correlating or 
extrapolating the performance of high-energy chemical rocket engines . To a large 
extent, most of the early work was concerned with defining simplified analytical 
procedures for assessing the point and ultimate extent of departure from chemical 
equilibrium in the gas flow through typical convergent-divergent nozzles (e.g., 

Refs. 1, 2, 3, 4 and 5)* A few investigators (Refs. 6, J, 8, 9 and 10) addressed 
themselves to the problem of constructing, by exact numerical means, the entire 
set of flow properties within one -dimensional nozzles although few results were 
obtained. 

Most recently, much more sophisticated numerical procedures have become 
available (e.g., Refs. 11 and 12) and much work of both an experimental and analyti- 
cal nature has been reported on the evaluation of nonequilibrium flow effects in 
a variety of engine/nozzle configurations (e.g., Refs. 13, l4, 15, l6 and IT). 
Indeed, exact numerical procedures are now available (Refs. 11, 12, 18 and 19) 
which solve simultaneously the chemical kinetic and gas dynamic equations pertinent 
to two-dimensional and axisymmetric flow fields, such as those encountered in the 
supersonic portions of typical expansion nozzles. In combination with these 
kinetic procedures, additional computational schemes have become available for 
treating the turbulent boundary layer (e.g. Ref. 20) and for handling the transonic 
flow problem (Ref. 21). It has thus proved feasible within the last few years to 
undertake detailed analysis and correlation of actual rocket engine test data and 
to predict "a priori" the delivered performance of propellant systems and engine 
configurations of future interest for space application (Ref. 22, 23, 24, 25, 26 
and 27 ) . The performance data which have been analyzed and correlated by this 
contractor are reported in Refs. 27, 28, 29, 30, 31 and 32. Generally, correla- 
tion of performance data has been achieved for the five propellant combinations 
presented in this handbook. Selected H2-F2 data available at low pressure and 
high area ratio (Ref. 32) did not correlate well with the large mass of data re- 
ported in Ref. 27, 28, and 29 ; the measured performance reported in Ref. 32 is 
higher than predicted delivered performance for similar operating conditions and 
nozzle geometries. 

Use of computational procedures such as those described above is not entirely 
straightforward, however, for several reasons. First, each procedure makes use 
of selective input data, restrictive assumptions and/or optional subroutines 
which result in significant variations in computed performance results; second, 
the performance losses predicted by procedures which individually handle such 
phenomena as nonequilibrium nozzle flow, viscous drag, and exit flow divergence 
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are not necessarily additive losses; and third, most procedures are lengthy and 
all require access to high-speed digital computation equipment. In order to 
develop effective methods to accommodate these problems, work has been performed 
by this Contractor (Ref. 25) and by the Performance Standarization Working Group 
of the Interagency Chemical Rocket Propulsion Group (ICRPG) (Ref. 33) which has 
resulted in defining "recommended" computer codes and both simplified and standard 
methodologies for predicting and correlating liquid rocket propellant engine per- 
formance. Remaining has been a more exact specification of the interaction between 
the combustor and nozzle components and development of an approximate method for 
estimating rapidly, without digital computation equipment, the deliverable per- 
formance of propellant combinations of current interest. It is for this latter 
purpose that the charts assembled herein have been prepared. 

Specifically, a procedure of verified accuracy has been developed which pre- 
sumes the various component losses in an expansion nozzle can be deducted from 
the theoretical one-dimensional, shifting-equilibrium performance of the nozzle, 
neglecting any coupling effects which may actually occur in the flow system. 

Using this procedure a series of performance estimation graphs has been prepared 
for five propellant combinations typical of current earth-storable, space-storable 
and cryogenic systems (viz. aerozine 5 O-N 2 O 4 , CH^-C^ (17.^)-F2(8 2 *6$).> B 2 Hg-GF 2 .> 
l^-Og and I^-Fg)* Included are kinetic performance and aerodynamic degradation 
charts which are sufficiently general to permit user-supplied nozzle contour data 
to be employed to estimate the nonequilibrium performance of virtually any nozzle. 
Additionally, sufficient information is provided to enable estimation of the 
effect on kinetic performance loss of variations in reaction rate constants re- 
lative to the values used in the handbook. 

Every attempt has been made to maintain performance prediction accuracy 
within ± 1 # of the performance predicted by the full-kinetic axisymmetric machine 
program (TDK) recommended by the ICRPG (Ref. 33)* 
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OUTLINE OF PERFORMANCE PREDICTION METHOD 


Work previously reported by the United Aircraft Research Laboratories (Refs. 
22 - 25 ) has shown that reliable prediction of deliverable rocket engine performance 
can be made provided that accurate and complete kinetic data are available to 
assess nonequilibrium flow effects in high area ratio nozzles and that appropriate 
accounting can be made of influences imposed by inefficient combustion and nozzle 
aerodynamics . This work combined with additional efforts reported in Refs . 1 6 , 

17 and 34 suggest that simplified procedures may be employed to approximate the 
nonequilibrium flow effects for common high-energy liquid rocket propellants and, 
further, that aerodynamic performance losses due to viscous drag and nozzle exit 
flow divergence are additive with nonequilibrium losses when the appropriate 
optional gas model is employed to evaluate the aerodynamic losses. Thus, exclu- 
sive of combustor-induced losses such as incomplete reaction or flow striation, 
all major sources of nozzle performance degradation can be treated separately. 

It is this fact which serves as the basis for the method described herein to 
estimate deliverable engine performance . 

In outline form, the prediction method involves a sequence of steps, each 
undertaken to evaluate one of the component performance losses relative to the 
theoretical, one -dimensional thermochemical performance of the propellant combina- 
tion under consideration. Charts are provided to permit interpolation of each 
component loss over wide ranges of operating conditions and for different nozzle 
geometries and engine sizes. The charts have been compiled recognizing two primary 
modes of utilization: a "performance" mode and a "contour selection" mode. Under 

the performance mode, the desired objective is to predict the net deliverable 
performance of a given propellant combination for a specified set of operating 
conditions within a prescribed engine/nozzle configuration. Under the contour 
selection mode (not necessarily an optimum contour) on the other hand, the de- 
sired objective is to arrive at a nozzle size or configuration which will produce 
a specified level of performance for a given propellant combination under a pre- 
scribed set of operation conditions. The logic for each mode is given below. 

The Performance Mode 

Specified under the performance mode are the following : 

1. The propellant combination 

2. Engine operating conditions (including 0 /f ratio and combustion 
chamber pressure, P ) 
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3. The nozzle configuration and size (including throat radius, r^; 
radius ratio, r c /r^.; exit area ratio, ( ) exit; and nozzle 
contour — i.e. conical nozzle with prescribed expansion angle, 
hell nozzle with prescribed truncation, or other generalized type 

with area ratio gradient, . specified as a function of 

dx 

local nozzle area ratio). 

To be determined from the appropriate performance charts and tables are : 

1. The theoretical, one -dimensional, shifting-equilibrium and frozen- 

flow performance ). 

&Jt Vac 

2. . The freezing-point location, denoting the location beyond. which 

recombination reactions are no longer possible kinetically. 

3. The one -dimensional, partial-equilibrium performance (determined 
by the freezing-point location), from which the kinetic loss can 
be inferred. 

4. Aerodynamic losses, including 

a. The divergence loss 

b. The viscous drag loss 

c. The transonic-flow loss (assumed to be 0.5 $ of the shifting- 
equilibrium performance for most common nozzle configurations) 

5. The net deliverable specific impulse performance, which is equal 
to the partial equilibrium performance (from 3 above) minus the 
divergence, drag and transonic aerodynamic losses (from 4 above). 

6. The net deliverable thrust (from 5 above and tabulated equili- 
brium throat mass flow rates). 

The Contour Selection Mode 

Specified under the Contour Selection mode are the following : 

1. The propellant combination and operating conditions. 

2. The desired thrust level 


4 



3. The general nozzle configuration (conical or hell nozzle) and 
overall expansion area ratio, (A/A m ^ n ) e xit * 

4. Geometrical constraints (r c /r^_ and expansion half -angle, 6 , for 
conical nozzles; design area ratio — i.e. the area ratio at which 
the exit flow becomes parallel — for bell nozzles). 

To be determined from the appropriate performance charts are : 


1 . 


The theoretical, 
flow performance 


one -dimens i onal , 
^ s Pyac ^ * 


shifting-equilibrium and frozen- 


2. The shifting-equilibrium specific mass flow rate, vr/A, through the 
•nozzle throat. 


3. The nominal size, r^, of the nozzle throat (from the specified 
thrust level and the specific mass flow rate, from 2 above). 

4. The freezing point location, denoting the location beyond which 
recombination reactions are no longer possible kinetically 0 

5» The one -dimensional, partial-equilibrium performance (determined by 
the freezing point location), from which the kinetic loss can be 
inferred . 

6. Aerodynamic losses including 

a. The divergence loss 

b. The viscous drag loss 

c. The transonic -flow loss (assumed to be 0.5 $ of the shifting- 
equilibrium performance for most common nozzle configurations) 

7* The net deliverable specific impulse performance, which is equal 
to the partial-equilibrium performance (from 5 above) minus the 
divergence, drag and transonic aerodynamic losses (from 6 above). 

8. The net deliverable thrust (from 7 above and the equilibrium 
throat mass flow rate from 2 above ) . 


As the net nozzle thrust is necessarily less than the specified nominal 
thrust, it is required that an adjustment be made to the nozzle throat size 
(calculated in 3 above) and the series of calculations from 4 through 8 repeated. 
The number of iterations necessary to achieve the desired net thrust may, of 
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course, "be reduced initially by increasing nozzle throat size slightly over the 
value calculated in step (3), 

Many variations, too numerous to specify, are possible on the above proce- 
dures when specific effects of various parameters, such as r c /r^. for conical 
nozzles or design area ratio for bell nozzles, are to be studied. To permit such 
studies, the charts have been compiled with as much generality as possible, and 
full exploitation of the data contained herein is left to the user. 
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DESCRIPTION OE TECHNIQUES 


In order to establish a more complete understanding and insight into the 
application and manipulation of the charts displayed in the various sections of 
this handbook, as well as to achieve a better realization of the limitations of 
the methods employed to calculate the various flow losses, a description of each 
of the component loss mechanisms is presented below. 

Kinetic Losses 

Kinetic losses are evaluated herein using the modified Bray criterion for 
multireaction systems (Ref. 17). In order to fully utilize the data presented 
in this handbook, a thorough understanding of this criterion is essential, al- 
though straightforward calculation of nonequilibrium performance for typical 
conical and bell nozzles does not require such familiarity. 

Description of Bray Sudden-Freezing Criterion 


The Bray criterion is an approximate procedure for predicting the point in a 
reacting nozzle flow where a reaction has departed significantly from equilibrium. 
This is accomplished by determining the point at which the forward rate of reac- 
tion becomes of the same order as the rate required to maintain equilibrium. This 
concept is illustrated by considering a nozzle expansion in which an arbitrary 
three-body recombination reaction takes place, 


2J + M ssSs- J p + M (1) 

kb 2 

where the change of J 2 concentration with time at constant density is found from 
phenomenological kinetics (Ref. 30 ), 



k f [J 2 1 [Ml 


\[jf rj] eq 2 } 
{ [J 2 ] " [Jgleq S 



( 2 ) 


In Eq. (2) the bracket, it indicates the actual instantaneous concentration of 
the indicated species, k f , denotes the forward rate constant, and the sub-script 
"eq" denotes the instantaneous equilibrium values. For near-equilibrium flow, it 
is clear that Y - Y e ^« Y e ^, so that 


a [j 2 ] \ 
at Ip 


<< Q y eq 


( 3 ) 


while for near-frozen flow Y Y, so that 

eq 7 


d [J 2 1 

at 


~ ay > > ay 


eq 


00 
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Therefore, from Eqs. (3) and (4) it can he deduced that an approximate criterion 
for the freezing point of the reaction can he defined according to 

“ °»eq ^ 0p eq (5) 

as previously stated, where the forward reaction rate expression, r^, is evaluated 
with very little loss in accuracy using equilibrium quantities up to the freezing 
point . 


The Bray criterion has been used successfully to analyze flows in which only 
one reaction is energetically and kinetically important. However, this method 
of solution considering only a single reaction may he valid only for a relatively 
few complex chemical systems. 

For the expansion of products of combustion in a nozzle, several reactions 
must he considered if an effective point of chemical freezing is to he found. 
Extension of Bray’s one-reaction criterion to multireaction one-equilibrium noz- 
zle flows is not necessarily straightforward, unless one reaction out of the many 
taking place has tied up most of the energy, in the combustion process, or unless 
all energetic reactions tend to leave equilibrium at nearly the same location in 
the nozzle. When this is not the case, the simple freezing criterion can at best 
bracket the actual nozzle flow parameters (by assuming that the fastest and slow- 
est of the many reactions control). 

Modified Bray Criterion for Multireaction Systems 


The concept of a composite freezing point is based on the assumption that an 
important rate-controlling species exists whose net rate of formation or deple- 
tion becomes very small as the reactions in which it participates become very 
slow, and that once this species has frozen at some point in the nozzle (which 
will be referred to as the "composite-reaction freezing point"), all remaining 
species can no longer react. The problem of finding a composite-reaction freezing 
point reduces to determining the kinetically, and thermodynamically, important 
species and reactions which must be studied. The importance of many of the exist- 
ing species can be evaluated on the basis of (l) their relative equilibrium com- 
position values, a mole fraction of 0.005 being a practical lower limit, (2) the 
energy release of the chemical reactions in which they participate, and (3) their 
role in effecting the chain-breaking steps in the overall reaction mechanism. 

The reaction rate that is necessary to keep an important species in equili- 
brium can easily be found knowing the equilibrium composition history as calcu- 
lated using the conservation equations, state relations, and the pertinent 
Guldberg-Waage laws of mass action (Refs. 31 and 32). The required time rate of 
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change of the concentration for component i due to reaction only, i.e., excluding 
expansion effects , is 


( a [Mj] \ 

r ow d(Xj/w) 1 

V at Jp,e q 

L dx J 


( 6 ) 


where = px j_/W, x^ is the mole fraction of species i, and W is the mean 

n 

molecular weight, £ • While equilibrium machine solutions employ the local 

i=l 

flow area to minimum flow area ratio (i.e., A/A^^ ) as an independent variable, 
the axial distance, X, is a more convenient variable since 


/ a [Mj] \ 

V dt /p,eq 



d(Xj/W) 


) J eq 


d( a/a min ) 

dx 


(7) 


becomes indeterminate at the minimum area (i.e., at A = A ^ n ) . Evaluation of 
Eq. (6) can proceed once the equilibrium properties and a nozzle geometry have 
been specified. Use of the equilibrium values for computing the required gradients 
will be accurate up to the point where the first important reaction freezes; but, 
since the species under consideration will usually appear in this reaction, this 
does not represent a severe limitation on the method when, one reaction is control- 
ling or all reactions freeze at the same location. 


In complete analogy with the simple one-reaction Bray criterion, it is 
necessary when treating a multireaction system to compare the required reaction 
rate, Eq. (6), with the total forward rate of reaction written in terms of all 
reactions that make substantial contributions to the overall rate of production 
or depletion of species. The forward reaction rate equation for the net rate of 
production of species is 


a 


[Mj] 

at 



( 8 ) 


where N is the total number of reactions, n the total number of species, k„ . the 

forward rate constant for the reaction, and a'. . and a". . the stoichiometric 

t.h +h 

coefficients of the i species in the j reaction for the reactants and products, 
respectively, when the reaction taking place in the exothermic direction is 
written as (Ref. 30 ) 


n 

£ 

i=i 


a y Mj 


k fi 


n 

£ 

i» I 


ay Mi 


(9) 
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When a single rate -controlling reaction exists, for example, the oxygen recombina- 
tion reaction in the expansion of air up to 6000 K, Eg.. (8) takes the simple form 


a [Mil 
at 



Mo 


( 10 ) 


where M- refers to the species in the reaction which, in addition to being 
thermodynamically important, has a large variation in composition through the 
nozzle to ensure accurate differentiation of the composition with respect to axial 
distance in Eg. (6). An effective freezing point for a multireaction nozzle flow 
can then be defined from Egs. (6) and (8) as 





n 

k f . v 
M 



( 11 ) 


In the situation described above where one dominant reaction exists, the right- 
hand side of Eg. (ll) is replaced by Eg. (10) so that this freezing criterion 
reduces essentially to the simple Bray criterion. In any event, it is clear that 
there can be as many freezing points as radicals or atoms to which Eg. (ll) is 
applied, but these points usually remain in close proximity to one another in 
practice, depending on the choice of reactions employed in the forward reaction 
rate eguation. 


The most successful results from this composite reaction effective freezing 
point analysis are achieved when only three-body reactions are considered in 
evaluating the right side of Eg. (ll), since bimolecular reactions in which one 
or more new radicals are formed for every radical consumed add very little energy 
to the flow after the three-body reactions have frozen.. 

Graphical solutions to Eg. (ll) are easily carried out for arbitrary flow 
systems, once the important reactions and species have been ascertained, by 
performing graphical differentiation of the eguilibrium concentration profiles as 
a function of axial distance for a preselected nozzle geometry, X = X(A/A ff|1n ) . 
Then, for a given set of forward rate constants and with all the eguilibrium pro- 
perties known as a function of area ratio, the left side of Eg. (ll) can be plot- 
ted against the right side over the range of area ratios . The point of crossing 
indicates where significant freezing has occurred, and is thus termed the freezing 
point . 


Method of Calculation of Freezing Point 


As was defined in Eg. (ll), the location of the freezing point is determined 
to be that position within a nozzle where the eguilibrium composition gradient and 
the kinetic composition gradient become egual. Although the kinetic composition 
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gradient is independent of the nozzle area ratio gradient, the equilibrium Com- 
position gradient is proportional to the nozzle area ratio gradient d(A/A min )/dx 
for a prescribed value of area ratio, A/A t|| . T 1 . Specifically, from Eq. (j ) 


Mi) 


dt 


EQUIL 


-{ 


d M. 


d (A/A 


d (A/A . ) 
' ' mm' 


EQUIL 


'Amin)) 
dx ) 


or alternatively, 

'd [m, 


dt 


EQUIL 


= ( d [mJ 

)d (a/a . ) 

( v ' min' 


EQUIL 




d ( A / A min)] 
^ (x/r t ) j 


( 12 ) 


(13) 


Therefore, for a selected propellant system, o/F ratio, and chamber pressure, the 
term 


d M. 


f d ( A / A min) 




EQUIL 


is independent of nozzle geometry. All of the geometric constraints particular 
to a given nozzle are in the second term of Eqs. (12) and (13). 

Use has been made of this fact in preparing generalized charts to permit 
determinations to be made of freezing points for each of the "five propellant 
combinations under consideration. In particular, both kinetic and "normalized" 
equilibrium gradients have been established as functions of local nozzle area 
ratio and plotted for several combustion chamber pressures and o/F ratios. Two 
normalized gradients are used; one which basically shows the dependence of the 
term 

{ a H v ) 

( d ( A / A min^ ) EQUIL 

on (A/A min ) by assuming 


d (A/A . ) 

, .. . . -.-33/ = 1 Ft" 1 

dx 


and the second which shows the gradient for a 15 deg conical nozzle with r / r^ = 1, 
r^ = 1.0 Ft. (0.3048m), The latter gradient provides a practical reference from 
which freezing points can easily be obtained for scaled 15 deg conical nozzles 


simply by multiplying the gradient displayed by the factor i_ . 

r t 


This former 
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gradient is the more general one and can he used to represent any given nozzle 
simply hy multiplying the gradient displayed hy the factor 


d ( A A mln ) 


dx 


' t 


Curves displaying the differential gradient. 


l d ( A / A min) 
d (x/r t ) 

3 (^-/Amin ) 
d "(x/r+) 


are provided for repre- 


sentative nozzle configurations in Figs. HD-1 through HD-4. Gradients are presented 
for conical nozzles of half angle from 5 to 30 deg and values of t c / t ^ up to 10 
in Fig. HD-1. Similarly, gradients for hell (perfect) nozzles are shown in Figs. 
HD-2 to HD-4. The gradient curves (Figs. HD-2 to HD-4) as well as corresponding 
nozzle profiles (Figs. HD-5 to HD-9) are presented for families of perfect nozzles 
with design- area ratios of from 100 to 1000, an r c /r^. = 1 and for 3 values of 
specific heat ratio, 1.25, 1.30 and 1.35. The data displayed allow the user to 
select a hell nozzle subject to any truncation criterion of his own choosing. 

The user is reminded that the equilibrium gradient to he employed to determine 
freezing point location is a function of throat radius and it is not the non- 
dimensional gradients which are illustrated in Figs. HD-1 through HD-4. 


As was indicated previously, the kinetic gradients displayed in the charts 
presented herein are independent of the nozzle area ratio gradients. However, 
the kinetic gradients are functions of the local equilibrium flow properties and 
compositions within the nozzle and the kinetic rate constants. For a simple 
recombination reaction such as that given hy Eq. (l), the appropriate kinetic 
gradient from Eq. (5) is 


3[Mj]\ 

dt //O.KIN 


{k f [Mj] 2 [M]}- 


(14) 


For the more general case involving multiple reactions as given hy Eq. ( 9 ), the 
appropriate kinetic gradient is 


d [m] \ 

d\ /^> KIN 



(15) 


where the term on the right-hand side of Eq. 15 is called the "composite" kinetic- 
gradient. 

In both Eqs. (l4) and ( 15 ), one recognizes that the kinetic gradient for 
each reaction is scaled directly hy the rate constant, k^, for that reaction. 

Thus, it becomes a relatively straightforward process to use the gradient charts 
provided herein to evaluate the sensitivity of the partial-equilibrium performance 
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to variations or changes in rate constants relative to those already employed. 

In the most simple case, the user may choose to work with the composite kinetic 
gradient only. Typically, in this instance, the kinetic gradient is scaled by a 
factor which is assumed to apply to the reaction rates for each individual reac- 
tion comprising the composite. In the more general case, however, factors can 
he applied to each individual reaction and a new composite gradient computed. To 
affect such changes the ratio of desired reaction rate to the rate used in the 
handbook can be determined from the reaction rate tables supplied herein for the 
specific propellant, and the individual kinetic gradient curves then scaled by 
these factors. The new composite kinetic gradient must be established by summing 
the contributions of the several reactions and the new freezing point then estab- 
lished. 

Graphical solution of Eg., (ll) to establish the freezing point results in 
defining an area ratio within the nozzle beyond which the energetic three -body 
recombination reactions cease to be important. The assumption is that no further 
chemical reaction takes place downstream of the freezing area ratio. This assump- 
tion may be validated by comparing performance results obtained on this basis 
with results of exact numerical calculations using a one -dimensional full-kinetic 
machine program such as those described in Refs. 12 and 33* In effect, such com- 
parisons calibrate the sudden-freezing criterion. For example, it was found in 
Refs . 22 and 25 that the effective freezing point for the H2-F2 propellant com- 
bination occurred at an area ratio where the composite kinetic gradient was equal 
to two times the equilibrium gradient. This delayed freezing apparently accounts 
for some continued energy release in the flow beyond the simple freezing point, 
which energy may be attributed to relatively fast two-body reactions (e.g., 

H2 + F — "-HF + H, in this particular case). This delayed freezing may be accom- 
modated by use of a transition factor in the sudden-freezing analysis which 
numerically equals the ratio of the composite-kinetic gradient to the equilibrium 
gradient required to achieve correlation between full-kinetic and sudden-freezing 
performance results . For the charts appearing in this handbook the following 
transition factors were employed. 

Propellant Combination Transition Factor 


Aerozine 5O-R2O4 
CH4-FI0X 


H 2 -°2 


h 2 -f 2 

b 2 h 6 -of 2 


1 

2 

1 

2 

2 
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Recent evidence (Ref. 2 6) indicates that a transition factor of 2 may he more 
appropriate for the CH4-FI0X and B 2 E£- 0 T ? 2 propellant combinations, and the user 
of this handbook is advised to adopt such a factor if the greatest precision is 
desired. The transition factors for the Aerozine 50 — and H2-O2 propellant 
combinations have been verified in Ref. 13. 

A tabulation indicating the location in this handbook where the kinetic and 
equilibrium charts can be found appears in the "Index to Figures and Tables". In 
this index, figures are numbered with prefix letters or roman numerals to refer 
to each specific section or propellant. For example, Fig. 1-13 to 1-30 provide 
the desired gradients for the H2-F2 combination; Fig. 11-13 to 11-30 provide 
those for the Hg-Og combination; Figs. III-13 to III-30 provide for the CH4-FI0X 
combination; Figs. IV-13 to IV-30 for the Aerozine 50-N2O4 combination. Figs. V-9 
to V-24 for 'the BgHg-C^ combination. 

Reaction Rate Data 


Of the five propellant combinations considered during the course of this 
analysis, viz. H2-0 2 , H 2 -F 2 , A5O-N2O4, CH4-FI0X (82.6$ F 2 ), and BgHg-CFg the 
CH4-FI0X and B2H5-OF2 are the most complex. The CH4-FI0X propellant combination 
involves at least twelve (12) gaseous species which can be related by twenty- 
three (23) reactions when all third bodies are considered as equally efficient 
(Ref. 38). The number of reactions is substantially higher if particular third- 
body efficiencies must be identified for the propellant combinations and also if 
condensed species, e.g. carbon were considered. The reaction mechanisms which 
describe the recombination process of the other propellants are subsets of the 
CH4-FI0X system, since the combustion species that are present for the H2-O2, 

H2-F2 and A5O-N2O4 (I2 excepted) are common to the species generated in GH4-FI0X 
combustion. The exception, nitrogen, which is present in A5O-N2O4, and possible 
nitrogen containing species are not involved in significant recombination reac- 
tions; however, the nitrogen and its compounds may be significant as third bodies. 

The reactions contributing significantly to the recovery of energy during 
recombination processes are usually limited to three-body reactions. Many other 
reactions, such as bimolecular reactions, occur in the full mechanism, but they 
can usually be neglected as secondary in establishing the nonequilibrium per- 
formance of the propellant combinations (Refs. l6, 22, 38 and 39)* Although the 
bimolecular reactions are important to the production of selected species for 
the three-body energetic reactions, their influence on propulsion performance 
is felt indirectly and only as long as the conditions are present in the nozzle 
for the three-body recombination reactions to occur at the minimum critical rate 
(freezing point) (Ref. 16). With the view of reducing the mechanisms to a toler- 
able number of reactions that can be treated utilizing the modified freezing- 
point analysis only significant three-body reactions can be considered for the 
complex combination of species possible for the CH^ fuel when combined with Fg-Og 



bearing oxidizers and the other fuel-oxidizer combinations. The mechanisms, 
therefore, have been limited to the principal recombination reactions : 

H + H + M - Hg + M 

H + F + M *- HF + M 

H + OH + M -HgO + M 

where M represents significant third bodies unique to the propellant combination. 
The selection of the particular three-body reactions has been based in part on 
analysis performed to correlate space-storable engine test data (Ref. 26). Also 
studies to establish the sensitivity of performance to reactions and rates in- 
clude the 'recent investigations by the ICRPG (Ref. 38) and RASA (Ref. 39). The 
studies corroborated the conclusions presented earlier in Refs. 13, l6 and 22. 

The BgHg-OFg propellant combination involves at least sixteen (l6) gaseous 
species and as many as fifty (50) reactions when all third bodies are considered 
as equally efficient (Ref. 38). As in the CH4-FI0X combination the number of 
reactions is substantially higher if particular third-body efficiencies must be 
identified. Analagous to CH^-Flox the important reactions are limited to the 
three body recombination reactions listed on page 19 with appropriate consideration 
given to 3rd bodies unique to the propellant combination namely, BF and 

BOF. The third body efficiencies for BF and BOF were estimated during a separate 
phase of the contract involving analysis and correlation of B2Hg-0F2 performance 
test data (Ref. 26). The other third bodies, namely H2, H, HF, and H2O are common 
to the propellant combinations . 

The recombination mechanism for each propellant combination has not been 
simplified to the point that only general third bodies, M, are considered. How- 
ever, the number of individual third bodies have been reduced to include only 
those of prime importance. It was assumed that the mole fraction concentration 
of a species must be at least 10“^ to be potentially important as a third body 
for the three-body recombination reactions. Species whose mole fraction concen- 
trations are less than 10" are considered only if they are prime reactants or 
products in the mechanism, e.g., OH in the water recombination reaction. 

This compromise approach between general and individual third bodies allows 
maximum flexibility in the kinetic performance charts so that individual refer- 
ence reaction rates used in this investigation can be changed when better data 
become available or can be changed to reflect differences in rates as reported in 
the literature. This is in contrast to the rates listed in the ICRPG compilation 
presented in Ref. 38, for which only a general third -body effect can be assessed. 
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The table below summarizes the reaction rate constants and the third-body 
efficiencies which can be combined to establish the reaction mechanism and rates 
appropriate for each of the propellant combinations , viz. Hq-Oq, S 2 -F 2 , A50-N204, 
CH 4 -FI 0 X and B 2 Hg- 0 F 2 - 

Summary of Elementary Reactions and Reaction 
Rate Constants 


Reaction 

Forward Rate: lbs -moles, 

ft -3 , sec, °R 

1 H + H + At • Ar 

kp = 4.62 x lO^T" 1 

2 H + F + Ar — nr HF + Ar 

k f = 1.155 x 10 15 T _1 

3 H + OH + Ar = 3 = H 2 0 + Ar 

k f = 7*85 x lO 15 ^ 1 


Third-Body Efficiencies Relative to Argon 

Reaction 


Third Body 

1 

2 

3 

H 

20 

1 

3 

h 2 

2.5 

2.5 

3 

HF 

2.5 

2.5 

3 

co 2 

2 

2 

10 

CO 

2 

2 

3 

h 2 o 

2.5 

2.5 

20 

h 2 

2.0 

— 

3 

BOF 

2.0 

2.0 

10 

BF 

2 

2 

3 


The H 2 -O 2 propellant combination mechanism includes reactions (l) and (3) 
with third bodies H, H 2 and H 2 O (a total of six reactions). The. H 2 -F 2 propellant 
combination mechanism includes reactions (l) and (2) with third bodies H, H^ and 
HF (a total of six reactions). 
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The A50-NgO^ propellant combination mechanism includes reactions (l) and (2) 
with third bodies E, HgO, CC> 2 > CO and N 2 (a total of 12 reactions). 

The CE^-Flox propellant combination mechanism includes reactions (l), (2) 
and ( 3 ) with third bodies E, Eg, BF, COg, CO and EgO (a total of 18 reactions). 

The B2Bg-0F2 propellant combination mechanism includes reactions (l), (2) and 
(3) with third bodies E, B 2 , HF, Eg 0 , BOF, and BF (a total of 18 reactions). 

Individual tables of the reactions and rate constants have been expanded for 
each propellant combination from the above summary and are included separately in 
the handbook compilation (c.f. Index to Figures and Tables). 

The rate constants for the Eg-Fg system are consistent with the literature 
data and have been selected by UARL as being representative of the most reliable 
data available at this time (Refs. 40 and 4l). 

The rate constants for the B 2 -C^ recombination reactions are also consistent with 
the literature data. The hydrogen recombination reaction: 

H + B + M H 2 + M 

with B and Eg as third bodies is common to both the Eg-Fg and Eg-Og propellant- 
combination. The rate constant for the water recombination reaction with argon as 
a third body: 


B + OB + E 2 0 + M 

has been established from analysis of the experiments of Schott and Bird (Ref. 42) 
along with the dependent hydrogen recombination rate for: 

B + H + M ===== B 2 + M 

as established in Ref. 40 (see e.g. Ref. 43). The effect of EgO as the third body 
on the E 2 O recombination is about 20 times that of argon as established by compari- 
son of the rate with argon (Ref. 43) and experimentally measured values in flame 
studies for which EgO was the significant third body (Refs. 44 and 45). 

The third-body efficiencies for E, Eg, BF, COg, CO, EgO, and Ng, BOF, and BF 
relative to argon for reactions not cited in literature have been estimated and, 
therefore, are not based on experimental verification. For example, the third- 
body efficiencies of CO 2 , CO, EgO and Eg for reaction (l) of the above table are 
estimated to be near that reported for Eg. The order and relative third -body 
efficiencies estimated for reaction ( 3 ) above are generally consistent with the 
values recommended in Ref. 46. 
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Equilibrium & Frozen Performance 


The calculation of the equilibrium, frozen, and partial-equilibrium flow 
properties and performance was performed with the machine program described in 
Refs. (36 and 37) using latest JANAF thermochemical data (Ref. 47). The Ref. 37 
computer program has been found to give results identical to those from the ICRPG 
reference one -dimensional equilibrium (ODE) program (Ref. 33). 

Estimates of the maximum performance that can be obtained from a specific 
propellant system can be obtained from the equilibrium performance curves presented 
for each of the propellant systems considered. An estimate of the maximum pos- 
sible kinetic performance loss can similarly be obtained by comparing the frozen- 
flow performance with the equilibrium performance for a particular propellant 
combination,* the difference is the maximum kinetic loss. 

The maximum thrust available from a particular propellant system with a 
nozzle of specified exit area ratio and throat radius is obtained from the equili- 
brium vacuum specific impulse and the appropriate equilibrium mass flow rate tabu- 
lated in this handbook (i.e. from Tables 1-2, II-2, III-2, IV-2, and V-2). The 
maximum thrust, T max , is then evaluated by: 

Tmax 7rr+. ( A )equilibrium ( Isp vac) equilibrium 

Additionally, the effective frozen specific heat ratio for a particular pro- 
pellant and o/F ratio is given in these tables for use with the bell nozzle design 
curves and divergence loss curves which require specification of specific heat ratio. 

Partial-Equilibrium Performance 


The procedures outlined in a preceding section of this handbook are utilized 
to obtain the freezing-point location for a given propellant system in a nozzle of 
specified geometry under a fixed set of assumed operating conditions. In order to 
use the freezing-point location to establish the corresponding nonequilibrium 
(i. e., partial-equilibrium) performance, it is necessary to perform thermochemical 
calculations which assume that shifting-equilibrium expansion occurs up to the 
freezing-point location, and that frozen-flow expansion occurs downstream of the 
freezing-point location. A large number of such calculations have been performed 
using the Ref. 37 (see also Ref. 33) computer program for each propellant combina- 
tion considered herein for broad ranges of o/F ratio, combustion chamber pressure 
and nozzle exit area ratio. The results of these calculations appear in graphical 
form to show the effect of freezing area ratio on partial-equilibrium specific 
impulse performance. These graphs are indexed in the "Index to Figures and Tables". 
Partial-equilibrium performance results for Hg-F^ H 2 -0 2 , CH^-Flox, A50-N 2 0^ and 
B 2 Hg- 0 F 2 , appear in Figs. 1-31 to -48, 11-31 to -48, IV- 31 to -48, and V-25 to 
V-40, respectively. Entering the appropriate curve with a freezing area ratio 
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and the nozzle exit area ratio, partial equilibrium specific impulse is obtained. 
It is to be emphasized that the specific impulse so defined must be reduced by 
the aerodynamic losses (to be described subsequently) in order to obtain the net 
deliverable performance. 


Aerodynamic Losses 

The primary aerodynamic losses which are encountered in expanding flows through 
convergent-divergent nozzles are three-fold: losses due to viscous dissipation 

along the nozzle walls., losses due to exit flow divergence and losses due to the 
nonuniform nature of the flow in the transonic region of the nozzle. 

Flow Divergence Loss 

The divergence losses presented in this handbook have been calculated by 
comparing analytical one -dimensional and axisymmetric thrust coefficients for the 
nozzle types, propellant combinations and conditions considered in this handbook. 

The computer programs used for these calculations include the two-dimensional 
equilibrium program (TDE) recommended by the ICRPG (Ref. 33) and developed by the 
United Aircraft Corporation (Ref. 48) and the one -dimensional equilibrium program 
described in Ref. 36 . 


Divergence losses are presented for families of perfect nozzles as a function 
of amount of nozzle truncation (percentage of the design area ratio) for three 
values of specific heat ratio (Figs. AL-1, -2 and -3). 


For conical nozzles, the divergence loss can be estimated for the larger exit 
area ratio nozzles — e.g. (A-/& m ± n ) exit = 20, as 


D.L. 


_ t 1-cos 9 ) ( T \ 

p s £Vae'frozen 


(IT) 


where 9 is the half -angle of the conical nozzle. Estimation of divergence loss 
using this expression is found to be within for exit area ratios greater than 
20; extension to lower values of area ratio may induce slightly greater divergence 
losses. Divergence loss is noted to be independent of nozzle scale. 


Additionally, divergence losses are presented (Fig. AL-4 and AL-5) for the 
15 deg conical and JO'fo bell nozzles at a series of pressure and thrust levels to 
allow immediate estimation of divergence loss for these commonly used nozzles. 
The TO'/o bell nozzles, as used here, refer to a truncated perfect nozzle 70 of 
the length of a 15 deg conical nozzle of the same exit area ratio. Thus, for a 
fixed specific heat ratio, a unique design area ratio is also implied when a 70 $ 
bell is discussed. The non-dimensional length of a 70$ bell nozzle of specified 
exit area ratio is given by the equation: 
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The factor 0.70 would be changed, of course, for the appropriate value for 
some other truncation criterion, say to (.60) for a 6of> bell, etc. 

X T 

The calculated coordinates of the exit point of a bell nozzle (Jn, _£.) are 

r t r t 

then plotted on the appropriate Figure: ND-5, -6, or -7 and the nozzle contour 

required of the design area ratio selected, usually by interpolation. Nozzle area 
ratio gradients are presented for most of these nozzles in Figure ED-2, -3 and -4 
and the appropriate nozzle area ratio gradient obtained for a particular nozzle by 
interpolation. For these later results divergence losses were evaluated using a 
chemically-frozen gas model and the corresponding option in the axisymmetric-f low 
machine program (Ref. 48). Experience has shown the effects of flow divergence 
and nozzle curvature can be accurately estimated by this flow model (Ref. 25). 

Boundary Layer Loss 


Frictional boundary layer losses have been evaluated for 15 deg conical and 
70$ bell nozzles with the modified Bartz boundary layer procedures (Ref. 20) 
using the ICRPG reference turbulent boundary layer (TBL) computer program (Refs. 33 
and 49). In using this program a frozen, axisymmetrical-flow gas model was assumed 
together with the adiabatic wall temperature option. Use of this gas model and 
wall temperature option has been shown to provide results consistent with results 
of rocket engine data correlations '(Ref* 25). 

Results obtained over wide ranges of conditions are presented in Figs . AL-6 
to AL-13. A sufficient number of pressure levels and thrust levels have been 
examined to allow interpolation of drag loss over the range of pressure and thrust 
values listed for each propellant in the figures. In contrast to divergence loss, 
boundary layer loss is known to be a function of nozzle scale. No effect of wall 
heat transfer has been included in the handbook. The user is referred to Ref. 33 
for a more comprehensive treatment of the boundary layer loss mechanism. 

Transonic Loss 


The flow characteristics in a sonic throat are such that the sonic line is 
more nearly parabolic than linear and that its position relative to the minimum 
cross-section is dependent upon the geometric parameters of the throat. It is 
evident that mass flow calculations based on the assumption of uniform axial flow 
across the minimum section will be somewhat in error. 
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The transonic loss accounts for nonuniform transonic flow and for frictional 
drag loss in the subsonic portion of the nozzle. For throat curvature to throat 
radius ratios, r c / r t> greater than 2 the procedures developed by Hall (Ref. 50) 
could be utilized. However, the use of such simplified procedures are not appli- 
cable for small values of throat approach (subsonic) r c /r.j. (e.g. r^r^ = 1). For 
such configurations the transonic solution developed at UAC (Ref. 21) is more 
applicable and was utilized in order to obtain a corrected mass flow rate, and 
impulse performance. On the basis of a large number of such calculations it has 
been found that the transonic loss amounts to a maximum of approximately 0.5 $ of 
the shifting equilibrium performance for reasonably configured nozzles (viz. 
r /r^ >0.5). Because of the very small magnitude of this loss, no charts have 
been prepared to permit its interpolation, and it is recommended that the user 
assign a conservative value of 0.5$ for all nozzles, propellant combinations and 
conditions' considered herein. 


Performance Estimation for Gaseous Propellants 

Tables and figures presented herein apply to liquid oxidizer-liquid fuel 
systems. Performance estimates for other than these liquid oxidizer-liquid fuel 
systems are determinable in the form of 4l s .p values to be used as corrections to 
the liquid-liquid estimates. A general equation, applicable to gaseous oxidizer- 
liquid fuel, liquid oxidizer-gaseous fuel and gaseous oxidizer-gaseous fuel 
systems is presented. 


Sufficient data are included in Table 1 to allow calculations of A I for the 

sp 

various propellants in their commonly used states. In a typical application, for 
example, cryogenic propellants are injected at their respective boiling points, 
hence specific heat data for the liquid phase of these components is not presented. 
In the A 50 -N 20 ^ system, on the other hand, Aerozine 50 is commonly injected at a 
"standard" temperature of 77°F (25°C) and as the boiling point of Aerozine 50 is 
considerably above this value, it is assumed that only liquid phase injection is 
of interest. Data for N 2 O 4 is included for both liquid and gaseous phase since 
the boiling point of this component is 70 C E, and the options of injection of the 
liquid at a temperature other than the boiling point or of gaseous injection are 
available to the user. 


An average specific heat has been presented for the propellants in Table 1; 
no great effect on 4l g was found to occur as a result of this approximation and 
it does afford considerable simplification in the calculation of 4I 0 „. 

sp 

The performance of a system with one or both components in a gaseous state 
may be estimated from the approximate equation: 


I 


sp 

9i 

i-g 


M 




(19) 
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where I 


S P„ 


Equilibrium vacuum specific impulse of the liquid-liquid systems; 


J = 778 £4*2- 3 ? . ; g = gravitational constant (32.174 ft/sec^) (in metric units, 
BTU 

j/g = 1 

kg-m 


and 


AH ^ H (r UE;i _ + (0 /F) (^ h oxidizer) 
I + 0/F 


( 20 ) 


where: ^ %uel = (Enthalpy of fuel at desired state) - (Enthalpy of liquid fuel 

at its boiling point) 

^ H oxidizer = (Enthalpy of oxidizer at desired state) - (Enthalpy of 
liquid oxidizer at its b.p.) 

0 /F = oxidizer-fuel weight ratio 

Since A I corrections are necessarily small relative to the total equili- 
brium performance value, calculations may be simplified by retaining only the 
first term of a Maclaurin series expansion of the full equation for approximate 
performance gain: 


Defining : 


one obtains : 


A ^sp Isp I sp. 
Q4 
t g 
g.g 


AI sp = U/2*)(^A.H)(— 


sp^, 


( 21 ) 


( 22 ) 


which can then be applied as a correction to the liquid-liquid performance obtained 
from the tables and figures presented in this handbook. 


Values of the heat of vaporization for each fuel and oxidizer used in the 
handbook are presented in Table 1 and are for the particular component in the 
gaseous state at its boiling point. Values of average specific heat are also 
presented to allow calculation of sensible enthalpy change of the propellant 
component if it is desired to inject the component into the combustion chamber 
at some temperature T other than its boiling point: 

dT 

•liquid ' 1 ' vapj j T 'Jg as 


AH 


B.P. j 


H 


injection 


,+r s. dT 

iJtc 


( 23 ) 
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where j = FUEL or OXIDIZER 


or, in terms of the average specific heat, C : 


AH J = C P) liquid ^ TBP j T i nject) + H vapj + c pj gas ^ T j Tep j ^ 


whereidEUn-n ; T-n-n and CL are obtained from Table 1. 

VAi-j Dbj Pj 

The resulting is then used in Eq. 20 to find the mean Si to use in 

Eq. 22 to establish the performance change due to the altered component. 

As an example, choose a Hg^-OgCg) system operating at an o/F = 6 with the 
gaseous oxidizer injected at 100°F above its boiling point. 


From Eq. (24) 


AH oxidizer = 91. 69 BTU/lb + 100 x C = 104 BTU/lb 

Jr 


From Eq. (22) 


Al S p 2 


m m 


_ 1 A8. 362(104)' 


I is then selected from the appropriate figure at the specified chamber 
sp ii 

pressure and nozzle exit area ratio, for example, from figure II-l, 

@ 100 psia, (A/A min ) ex . t = 100, 

I = 465.2 sec (4559*0 N’sec/kg) 
sp f £ 

Al 0 ^ = 5.4 sec (52.9 N*sec/kg) 


Net performance for the HgO J-OgCg) is then: 

I st) = I + Al st) = 465.2 + 5.4 = 470.6 sec (4611.9 N* sec/kg) 

y Hg 
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Table 1 


Thermodynamic 

Properties 

for Performance 

Estimation of 

Gaseous 

Propellants 


^^vaporization 

T b oiling 

Cp 

Component 

Btu/lb 

Kcal/mole 

°F 

°C 

Btu/lb *°R or gSSL. 






^g. 0 

°2 

91.7 

1.63 

-297.4 

- 183.0 

0.22 

*2 

191.8 

0.2148 

-423.2 

-252.9 

3.50 

F 2 

73-9 

1.56 

-306.7 

- 188.1 

0.19 


219.4 

I .96 

-242.5 

-152.5 

0.52 

Flox 

77.0 

1.572 

-305.0 

-187.3 

0.18 

N 2°4 

178 . 2 * 

9.11 

70.0 

21.15 

O. 37 ** 






0.65 

A50 

425.8* 

9.54 

190.0 

87.8 

0 . 70 ** 

of 2 

88.67 

2.66 

-229.5 

-145.3 

0.17 

b 2 h 6 

221.2 

3.41 

-134.5 

- 92.5 

0.36 


* at 77°F (25°C) 

** Liquid 

Summary of Relevant Equations for Performance Calculations 


max 


IT 


1 1 V A /equilibrium 


UsD..,J 


vac 


equilibrium 


(16) 


Conical nozzle divergence loss 


•cos 8 


( Is Pvac ) 


frozen 


(IT) 
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Xe 


rt 


= 0.70 




0.26795 



(18) 


sp 
9 'i 

g-g 


T- ■ AH + x «t t 


(19) 


AH = 


AH Fuel + 0/F (^H oxidizer ) 

I + 0/F , 


( 20 ) 


AX_ = Icn * 1 


_ sp S P 


g-£ 
t g 
g-g 


s p, 




( 21 ) 


Xsp^i^A-H) 


-SP, 


jf-i- 


( 22 ) 


/» T b pi /* T j 

AH i = X C Pi,; quid dT + X- 0 ?*- 

'inject bpj 


(23) 


ah, = c r 


T or , - T 


i " W PJ liquid ^' BP j ' inject ) + H vapj ' "Plgas 


+ C, 


( T j T BP; ) (24) 


Transonic LOSS - 0.005 (^sp V ac ^equilibrium 


(25) 
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( Isp vac) net = ( Isp voc) portlol 2 ( AI sp)o«od,nGmL 

equilibrium 


T net (' I ' s Pvac ^net * ( A Equilibrium 


■A Equilibrium * ^ ( Xsp vac Equilibrium 


im 

' A ''equilibrium 


* ir * (I, n ) . 

v S P vac 'net 
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SAMPLE (SIMULATIONS 


Examples of the calculation procedure are presented to enable the user to 
familiarize himself with the mechanics of estimation of nozzle performance. 

Example 1 . Performance Estimation Of Specified Conical Nozzle 

Using the handbook in the performance mode, estimate the net performance of 
the following 15 deg conical nozzle operating as specified: 

(l) Nozzle: 15 deg conical 



(1) Equilibrium performance (Fig. I-l); frozen performance (Fig. 1-2); and 
maximum thrust (Eq. 16) and Table 1-2 

(2) Freezing point (Fig. I-l4) ; partial-equilibrium performance (Fig. 1-32) 

(3) Aerodynamic losses: Flow divergence loss (Eq. 17); drag loss (Fig. AL-6) 

and transonic loss (js'fo of equilibrium performance) (Eq. 25) 

(4) Net performance (partial- equilbrium I g p- £ (aerodynamic losses)) (Eq. 26). 

Actual determination of values for the above quantities proceeds as follows: 

(l) For this H2-F2 system operating at 100 psia and an 0/F = 12, and a 

nozzle of (a/A . ) . , = 100, the equilibrium performance is 488 sec 

K ' min' exit 

(4782 -* SeC ) obtained from Fig. I-l. Frozen performance is likewise 
kg 

obtained from Fig. 1-2 for this system and is found to be 415-5 sec 

(4072 N' sec ) . The partial equilibrium performance must lie between 
kg 

these two values. The maximum thrust obtainable with this nozzle 
can be found by using the equilibrium mass flow rate in Table 1-2, 
the equilibrium vacuum specific impulse and the nozzle throat radius 
in Eq. l6: 
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■*max 

7Tr+ 2 * T * (l S r) ) 

■ A- s £\rac ec L 

or 



T 

max 

3.1416 (0.3333) 2 (58.3) (488) 

or 

9931 lb (4.417 x 10%) 


%ax = 


(2) In Fig. I-l4, dividing the equilibrium gradient for the 15 deg conical 
nozzle with r^ = 1.0 ft (0.3048 m) by the specified throat radius 
r-fc = 0.3333 ft (0.1016 m) gives the new equilibrium curve for the 
smaller nozzle. The intersection of this new equilibrium curve 
with the composite kinetic curve establishes the freezing point 
which is found to occur at (A/ A^n ) = 2.10. 

An alternative procedure to establish the freezing point location is 
to utilize the second equilibrium gradient curve presented in Fig. I-l4 
along with the appropriate data from Fig. ND-1. 


For several values of area ratio, values of area ratio gradient are picked 
from Fig. MD-l. The selected points are given in columns 1 and 2 in table below. 
These area ratio gradient values are then divided by the throat radius of the 
desired nozzle and are tabulated in column 3 of the table. From Fig. I-l4 one 
obtains the kinetic gradient for a nozzle of unity area ratio gradient for each 
previously selected value of area ratio (col. l) and these are given in column 4. 
The net equilibrium H-atom gradient is then calculated by the following equation 
and listed in column 


or 


a [hi 
at 


1 


/ aH 

/ \ 

eq r t 

*• d^x /r t j 

1 at 7 eq. 

’ d r /A MiN/ ’ 


L dx ■ 'J 


aH 

at 


eq 


(Column 3) * (Column 4 ) 


and is calculated for each value of A/A^.^ in column 1. 

This resulting data is then plotted on Fig. I-l4 (an overlay should be used) 
and the intersection of the resulting curve with the "composite kinetic (transi- 
tion factor = 2)" curve in this figure establishes the freezing point. 
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4 


5 


l 


2 


3 


A/ A min 

f d( A//A MIN } 
V djx/r t j , 

1.02 

0.18 

1.05 

o. 4 o 

*1.08 

0.555 

1.20 

0.585 

1.50 

0.660 

2.00 

0.760 

2.50 

0.85 



0.54 

1.20 

1.665 

1.755 

1.980 

2.280 

2.55 


ec T d A / A MlNj . , 

( iH ) 
^ at > 

0.90 

0.486 

0.60 

0.720 

0.50 

0.833 

0.30 

0.526 

0.13 

0.257 

0.075 

0.171 

0.047 

0.120 


eq 


This procedure is necessary when other than 15 deg conical nozzles 
are examined. The values of area ratio gradient may be obtained 
for conical nozzles of various r c /r.j- and half angle specification 
from Fig. ED-1, and perfect nozzles of specified design area ratios 
from Figs. ED-2, -3, -4. For the conical nozzle-circular arc 
combinations shown in Fig. HD-1, the area ratio gradient is 
obtained by following the desired r c /r^- curve from the throat to 
the intersection of this curve with the conical section curve of 
the specified half-angle and then following this curve. The 
intersection of the two curves represents the point at which the 
circular-arc and conical section are tangent and the contour is 
conical downstream of this point. Using either procedure, once 
the freezing point is established, the partial equilibrium perfor- 
mance is then obtained from Fig. 1-32. For this case with 
(A/A- n )f r = 2.10, the nonequilibrium performance is 461.7 sec 

(4525 li?- ec ) . 
kg 


(3) For a 15 deg conical nozzle, divergence loss is calculated from 
Eq. 17 (or found from Fig. AL-4) . 


D.L. 


1-COS Q 
2 

1-0.96593 

2 


s ^vac frozen 

*415.5 


D.L. = 7.1 sec (69 — — ) 

kg 

By interpolation, the drag loss is estimated from Fig. AL-6 as (1.2$)* 

N • SGC 

^spv-ac Wen < H S- ^ “ ^*^-5 = 5.0 sec (48.9 -£-) 

for this nozzle. Transonic loss is M, * (l„^ ) . 

sp V ac equilibrium 


*Tangent point and transition from circular-arc throat to conical skirt 
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(4) 


. . . , N-sec. 

(Fig. I-l) and is 2.4 sec (23.9 ). 

kg 

Net performance is then the partial equilibrium performance (Fig. 1-32) 
less the component losses (Eq. 26): 


^sPvac^net ^ 1 s Pvac ^ P ar t i al- e q . ^ ^ ^sp^ aerodynamic loss 


= 461.7- (7.1 + 5 + 2.4) 


thus 


(5) 


<Isp vac4et 


= 447.2 sec 


(4383 


N«seC ' 

kg 


Net thrust is the net vacuum specific impulse times the equilibrium 
mass flow at the throat (Table II ) (Eq. 10) 


T net - < Wnet * <?> 


* A t 
A eq t 


T = 447.2 sec * 58.3 * tt * (3333) 2 = 9101/lb 

net sec 


(4.048 x lO^r) 


Performance results and the appropriate figures and equations are summarized 
below:* 


Equilibrium Performance 

488 sec (4782 N ’ sec ) 

kg 

Figure I-l 

Frozen Performance 

415.5 sec (4072 N.sec) 

kg 

Figure 1-2 

Freezing Area Ratio 

2.10 

Figure I-l4 

Kinetic Performance 

461.7 sec (4525 ^g|- c ) 

Figure 1-32 

Equilibrium Mass Flow 

58.3 lb/sec/ft 2 (284.7 k f J 

sec 2 m 2 

Table 1-2 

Maximum Thrust Possible 

9931 lb (4.417 x 10^ l^sec) 

kg 

Equation 16 

Aerodynamic Drag 

5.0 sec (48.9 N - sec ) 

kg 

Figure AL-6 

Divergence Loss 

7.1 sec (69 IjL|®£) 

Equation 17 

Transonic Loss 

2.4 sec (23.9 ^l 22 ) 

Equation 25 

Net Delivered Performance 

447.2 sec (4383 N ’ sec ) 
kg 

Equation 26 

Net Delivered Thrust 

9101 lb (4.048 x 10^ N-sec) 

kg 

Equation 27 


*A summary of relevant equations may be found on p. 24. 
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Example 2. Performance Estimation Of Specified Bell Nozzle 


Using the handbook in the performance mode, estimate the net performance 
of the following hell nozzle operating as specified: 

(1) Nozzle; Bell, (A/A m ^ n ) ( ^ es ^g n = 300 
r c / r t ■ 1 -° 

r t = 0.3333 ft. (0.10l6m) 

(A/A . ) . = 100 

' min exit 

(2) Propellant (i) 

and 

Conditions: 0/F = 12 

P c = 100 psia (6.895 x lO? N/m 2 ) 

Quantities to be determined are: 

(1) Equilibrium performance (Fig. 1-1), frozen performance (Fig. 1-2) 
and maximum thrust (Eq. l6) and Table 1-2. 

(2) Freezing point (Fig. I-l4), partial equilibrium performance (Fig. 1-32) 

(3) Aerodynamic losses: Flow divergence loss (Fig. AL-3)> drag loss 

(Fig. AL-8) and transonic loss ( 2 $ of equilibrium performance) 

(Eq. 25). 

(4) Net Performance (partial equilibrium I S p -• ^/aerodynamic losses)) 

(Eq. 26). 

Determination of values for the above quantities proceeds as follows: 

(l) For this H0-F2 system operating at 100 psia and an 0/F = 12 and a 
nozzle of (A / A m j_ n ) exit = 100, the equilibrium performance is 488 sec 

(4782 , ?. c ) obtained from Fig. 1-1. Frozen performance is obtained 

kg 

from Fig. 1-2 for this system and is found to be 415.5 sec (4072 — ° sec ) 

kg 

The partial equilibrium performance must lie between these two values. 
The maximum thrust obtainable with this nozzle can be found by using 
the equilibrium mass flow rate in Table II, the equilibrium vacuum 
specific impulse and the nozzle throat radius in Eq. l6: 
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or 


T v = n * r t 2 * 2 * (Ig-D ) o „ 
max ^ a Bp vac e< l 

T^ = 3.1^16 (0.3333) 2 * (58.3) * (488) 


or 


t = 9931 16 (4.417 x loSsr) 

max 

(2) The appropriate specific heat ratio is found, from Table 1-2 to he 
1.34, thus, the perfect nozzle contour from Pig. HD- 7 ( Y = 1.35) 
and the corresponding nozzle gradient from Pig. HD-4 ( y = 1.35) 
are utilized. Interpolation between the values of area ratio 
gradient for Y =1.30 (Pig. HD-3) and Y - 1.35 (Fig. HD-4) has 
been avoided since the actual specific heat ratio and the Pig. 

HD- 4 value differ by a small amount for the sake of clarity in 
this example. 

In order to establish the freezing point, the alternative procedure 
described in step (2), example 1 must be followed. The area ratio 
gradient values obtained from Pig. HD-4, the normalized equilibrium 
gradient from Pig. I-l4 and the calculated equilibrium gradient shown 
below: 



1.2 

1.07 

3.21 

0.30 


.963 

1.4 

1.32 

3.96 

0.17 


.673 

1.6 

1.44 

4.32 

0.12 


.518 

1.8 

1.55 

4.65 

0.091 


.423 

2.0 

1.65 

4.95 

0.075 


.371 

2.5 

1.84 

5.52 

0.047 


.259 


where 

0 rt 

1 r d(A/A M[N ^ 

(M) 

ci(a/A ) 1 



a t 

eq r t *> d/x/r^W 

V at 'e q; 

dx J 



The values of equilibrium gradient tabulated above are then plotted 
on Fig. I-l4 and the intersection of the equilibrium and composite 
kinetic curve employing a transition factor of 2 curves establishes 
the freezing point at (A/ Am-in ) = 1.75. 

The partial equilibrium performance is then obtained from Fig. 1-32. 
For the present example and the freezing area ratio of 1.75? the 


34 



IT * sec 

partial equilibrium performance is 457 sec (4479 — ) . 

(3) By interpolation from Fig. AL-3, the divergence loss is found to 
be 1.3$ times (l SPvac ) frozen (Mg- 1-2): 

D.L. = .013 • 415-5 

D.L. = 5.4 sec (SP.Q ff_‘ sec ) 

kg 

Interpolation on Fig. AL-8 gives an estimate of drag loss for this 
nozzle as 1 ' % * ( I sp vac ) frozen ( Fi S* I " 2 ) = - 01 * ^5-5 = 4.2 sec 

(4l - sec ) . As the 70$ bell nozzles had about the same truncation 
kg 

fraction as the present example, the 70$ bell curves are used for 
drag loss estimation without appreciable error. Transonic loss is 

¥!o (I m ) (Fig. i-l) = .005 * 488 = 2.4 sec (23.5 jL’.ggS). 
ap vac eq kg 

(4) Net performance is then the partial equilibrium performance less 
the component losses (Eq. 26) 

^spvac^net " ^ I sp vac ^ partial eq. ~ ^ (^^sp) aerodynamic loss 
= 457 - (5.4 + 4.2 + 2.4) 

< Wnet - ^ ~ (4361 *«) 

Performance results and the appropriate figures and equations are summarized 
below: 


Equilibrium Performance 

488 sec (4782 N * sec ) 

kg 

Figure I-l 

Frozen Performance 

415-5 sec (4072 N>sec ) 

kg 

Figure 1-2 

Freezing Area Ratio 

1.75 

Figure I-l4 

Kinetic Performance 

, __ ,, , N*sec. 

457 sec (4479 ) 

kg 

Figure 1-32 

Equilibrium Mass Flow 

58.3 St ft 2 ( 284 -7 ~ u) 

sec rt sec m ^ 

Table 1-2 

Maximum Possible Thrust 

9931 lb (4.417 x 10%) 

Equation 16 

Aerodynamic Drag 

4.2 sec (4l N * sec ) 

kg 

Figure AL-8 

Divergence Loss 

5.4 sec (52.9 N * sec ) 

kg 

Figure AL-3 
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Transonic Loss 2.4 sec (23.5 - ,s - c ) Equation 25 

kg 

Net Delivered Performance 445 sec (4361 se .. c ) Equation 26 

kg 

Net Delivered Thrust 9056. 16 (4.028 x lo\f) Equation 27 

Example 3. Selection of Conical Nozzle of Specified Thrust 

Using the contour selection mode, size a 15 deg conical nozzle of 
specified exit area ratio and throat radius ratio to deliver 1000 lb (4448N) 
thrust utilizing the specified propellant and operating conditions. 

(1) Nozzle: 15° conical 

(V Amin) exit - ^ 
r c /r t =1.0 

Deliver 1000 lb (4448N) net thrust 

(2) Propellant ^(.O-C^-O 

and 

Conditions: 0/F = 6 

P c = 200 psia (1.379 x 10 6 N/m 2 ) 

In order to establish the required nozzle throat radius the following 
quantities must be determined: 

(1) Equilibrium performance (Fig. II-3) and frozen performance 
(Fig. II-4) . 

(2) Equilibrium mass flow at the throat (Table II-2) and the approximate 
throat size: 

I 

(Eq. 28) 


(3) Freezing point (Fig. II- 17) and nonequilibrium performance (Fig. 11-35) « 

(4) Aerodynamic Losses: Flow divergence loss (Eq. 17), drag loss 

(Fig. AL-7) and transonic loss (J$ of equilibrium performance) 

(Eq. 25). 

(5) Since the nozzle throat size required to deliver a specified thrust 
is dependent on the aerodynamic and kinetic losses which themselves 
are functions of nozzle throat size, the exact nozzle radius must 
be established through an iterative process. 
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Evaluation of the above quantities follows: 

(1) For an H2-O2 propellant at an 0/F = 6 and P c = 200 psia and a nozzle 
of exit area ratio 60, the equilibrium performance from (Fig. II-3) 

is found to be 457.5 sec (4484 — sec ) . Similarly, the frozen 

kg 

performance from Fig. II-4 is 425 sec (4l65 sec -) . 

kg 

(2) Equilibrium mass flow at the throat is found from Table II-2 to 
be 124.2 lb/sec/ft 2 (606.5 — ) the approximate throat radius 
is then (Eq. 21) 


r t = 


(3) Once the throat radius 
be established. In Fig. 11-17, the equilibrium gradient for the 
15 deg conical nozzle with rt = 1.0 ft (0.3048m) is dividied by 
the throat radius determined above, r^ = 0,075 ft (0.0229m). This 
new equilibrium curve is then plotted on .Fig. 11-17 and its inter- 
section with the composite kinetic curve establishes the freezing 
point which is found to occur at (A/A m j_ n ) = 1.24. The nonequilibrium 
performance is then obtained from Fig. 11-35 and is 44l sec 

(4322 5^22). 

kg 

(4) For a 15 deg conical nozzle divergence loss is calculated from Eq. 

17: 


ecu ui- 


T 


2 . 7T • (I 


SP. 


vac 


eq. 


1000 




, (124.2) (3.l4l6) (457. 5)i 


r, = 0.075 ft (0.0229m) 


has been calculated the freezing point may 


D. L. 


D. L. 


cos 0 
2 

-0. 96593 


( I spyac) frozen (klg. 

) * 425 


( '~7~\ 


N-sec, 


II-4) 


Drag loss is obtained by interpolation from Fig. AL-7 and is 

1.4 %* (Io-n ) (Fig- II-4) = (0.014) (425) = 6.0 sec (58.3 ' - ) . 

■^vac frozen kg 
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Transonic loss is 0.005 * (l„_ ) . (Fig. II-3) = (0.005) (457-5) = 

s Pvac eq 

, , N*sec N 

2.3 sec (22.4 k "-') • 

(5) Net performance is then the partial- equilibrium performance less 
the component losses (Eq. 26) 


(I 


sp 


) = (I 


vac net 


s P V ac ^ partial eq. 


“ L(AI) 


sp aerodynamic 
loss 


= 441 - (7-2 + 6.0 + 2.3) 


Thus 

ttsivJ net - te5 ' 5 < 4170 ^ 

The delivered thrust is thus, from Eq. 20, 

^net = ^spya^net ^ eq n r t 2 

= (425-5)(124.2)(3.i4i 6)(0.075) 2 
= 933-9 lbs (4154N) 


As this delivered thrust is less than the specified thrust, throat 
radius should be increased slightly and the calculation repeated 
for this increased nozzle size and iteration continued until the 
required delivered thrust is obtained. A good approximation to 
the required throat radius may be obtained by using the calculated 
delivered specific impulse instead of the equilibrium value initially 
used in the radius equation (Eq. 22) : 

T \i 

^ n (-*-sp vac )net J 

For the present example, this becomes 

r = ( 1000 

* \(124.2)(3.i4i6)(425. 

r = 0.0776 ft (0.0237m) 

*c 

and is accepted as a final value for r^ since no discernible change 
in partial equilibrium performance or in the aerodynamic losses was 
found. 
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Performance results and the appropriate figures and equations are summarized 
below for the first iteration. 


Equilibrium Performance 
Frozen Performance 
Freezing Area Ratio 


, „ , R*sec x 

457-5 sec (4490 — ) 

kg 

425 sec (4l60 ) 

kg 


1.24 


Figure II-3 
Figure II- 4 
Figure 11-17 


Kinetic Performance 
Equilibrium Mass Flow 
Aerodynamic Drag 
Divergence Loss 
Transonic Loss 
Ret Delivered Performance 
Ret Delivered Thrust 


439 sec (4303 

kg 

124.2 — -- " - U (606.5 
sec rt 


_k£. 


sec m 1 


2 ) 


„ ,_ n „ R«sec N 

6.0 sec (58.3 1 ) 

7.2 sec (71.0 ■ % ! — -) 

kg 

2.3 sec (22.4 --, sec ) 

kg 

423.5 sec (4150 

kg 


Figure 11-35 
Table II- 2 
Figure AL-7 
Equation 17 
Equation 25 
Equation 26 


929.5 lbs (4134R) 


Equation 27 


Example 4. Establishing The Effect Of Reaction Rate Change 


Utilizing the handbook in the performance mode, establish the effect 
on kinetic performance of a 4 fold increase in the reaction rate constants 
for the specified propellant, operating conditions and nozzle: 


(1) Rozzle: 15 deg cone 

(A/Amin)exit = 100 
r t = 1.0 ft (0.3048m) 

r c/ r t = 1 

(2) Propellant 1^(1 )-Fp(f ) 

and 

Conditions: O/F = 8 

P c = 100 psia (6.895 x lQ^ R/m?) 

(3) Reaction Rates: Reaction rates uniformly adjusted to 4 times the 

reference rates used for the handbook preparation (Table I-l) . 
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Quantities to be determined: 

(1) Freezing point (Fig. 1-13) with standard rates (Table I-l) and 
corresponding nonequilibrium performance (Fig. 1-31). 

(2) Freezing point (Fig. 1-13) with new reaction rates (each rate in 
Table I-l times 4) and corresponding nonequilibrium performance 
(Fig. 1-31) • 

(3) Net Performance change due to reaction rate increase. 

Calculations of the above quantities follows: 

(1) The equilibrium gradient curve for the 15 deg conical nozzle, 
r^ = 1 ft (0.3048m) is shown in Fig. 1-13* The freezing point 
located at the intersection of the composite kinetic and 
equilibrium curves is at A/A ^^ = 1.90. The corresponding partial 
equilibrium performance is then, from Fig. 1-31, 471 sec 

{k6l6 ^ 22 ) . 

(2) Since all the component reactions have been multiplied by a factor 
of 4, the sum of the component reactions has also multiplied by 

4 and the composite kinetic curve is increased by this factor, and 
replotted on Fig. 1-13, establishing the new freezing point at 
A/Amin = 2.5* The corresponding nonequilibrium performance is seen 

from Fig. 1-31 to be 476 sec (4665 ' N * Sec ) . 

kg 

(3) A performance gain of 5 sec (49 1 * sec ) is thus realized by an increase 

kg 

in the reaction rates of a factor of 4. 

Performance results and the appropriate figures and equations are summarized below: 


Original Freezing-Area Ratio 

1.90 


Figure 

1-13 

Kinetic Performance 

471 sec 

/, N‘seC\ 

(i^e ^ ) 

Figure 

1-31 

Revised Freezing- Area Ratio 

2.5 


Figure 

1-13 

Revised Kinetic Performance 

476 sec 

,, „„ N-sec v 
665 ^ > 

Figure 

1-31 

Net Performance Change 

+ 5 sec 

(+4q N-secj 

kg 
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LIST OF SYMBOLS 


a 

A 

•^min 

Is Pyac 

J 

k f 

k b 

M 

o/f 

P c 

r 


r t 

t 

T 

V 

W 

w 

W/A 

x 


x. 


y 


Forward rate parameter (see Eq. 2) 

Area 

Minimum area 

Vacuum Specific Impulse 

General Molecular Species (see Eq. l) 

Reaction rate constant in forward direction 
Reaction rate constant in backward direction 
General molecular species or general 1 third body 
Oxidizer-fuel ratio 
Chamber pressure 
Radius 

Throat radius of curvature 

Throat radius 

Time 

Thrust 

Velocity in axial direction 
Molecular weight 
Mean molecular weight 
Specific mass flow rate 
Axial distance 
Mole fraction 

Concentration ratio defined by Eq. (2) 


4l 





a 


ij 


y 

o 

p 


Stoichiometric coefficients of the ^ species in the reaction 
for the reactants and products 

Specific heat ratio 

Conical (divergence) half- angle 

Density 


Subscripts 

eq Equilibrium conditions evaluated on the basis of infinite reaction 

rates. 

fr Frozen condition or freezing area ratio 

kin Conditions evaluated on the basis of finite reaction rates 

* Sonic conditions 

i, k Refers 

0 Refers 


to species i or k 
to reaction j 
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NOZZLE AREA RATIO GRADIENT, d (A/Ami I* 

d(x/r.) 
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FIG. ND-1 


COMPARISON OF NOZZLE AREA RATIO GRADIENTS FOR TYPICAL CIRCULAR 

ARC AND CONICAL CONTOURS 

r c /r f = RADIUS OF CURVATURE/THROAT RADIUS FOR ARC CONTOUR 
6 = HALF-ANGLE OF CONICAL CONTOUR 































NOZZLE AREA RATIO GRADIENT 
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FIG. ND-2 


NOZZLE AREA RATIO GRADIENTS OF SELECTED PERFECT NOZZLE CONTOURS 


SPECIFIC HEAT RATIO = 1.25 
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FIG. ND-3 


NOZZLE AREA RATIO GRADIENTS OF SELECTED PERFECT NOZZLE CONTOURS 


SPECIFIC HEAT RATIO = 1.30 

r / r = 1.0 
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AXISYMMETRIC PERFECT NOZZLE CONTOURS 

SPECIFIC HEAT RATIO = 1.25 
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FIG. NO-5 
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FIG. ND-6 
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TABLE 1-2 

EQUILIBRIUM AND FROZEN FLOW PROPERTIES AT 
NOZZLE THROAT 
H 2 (i) - Fg(i) 



P c -psia 


OXIDIZER - 

FUEL WEIGHT RATIO 



4.0 

6.0 

8.0 

10.0 

12.0 

l4.o 

16.0 

Equilibrium 

100 

55.4 

55.8 

56,9 

57.7 

58.3 

59.0 

59.8 

Mass 

200 

110.7 

110; 9 

112.8 

CO 

1 — l 
1 — i 

115.5 

116.7 

118.3 

Flow - 

300 

165.9 

165.9 

168.4 

170.6 

172.3 

174.0 

176.3 

w/a 

500 

276.4 

275.5 

279.0 

282.4 

285.1 

287.9 

291.5 

Lbm 

750 

414.6 

412.1 

4 16. 8 

421.5 

425.4 

429.4 

434.6 

Ft 2 . Sec 

1000 

552.7 

548.6 

554.2 

560.0 

565.0 

570.2 

576.9 

Frozen 

100 

56.4 

57.8 

59.3 

60.3 

61.3 

62.2 

63.3 

Mass' 

200 

112.4 

114.5 

117.1 

119.1- 

120.9 

122.7 

124.8 

Flow - 

300 

1 

168.1 

170.8 

174.4 

177 A 

180.0 

182.6 

185.6 

w/a 

500 

279.5 

282.8 

288.3 

293.0 

297.I 

301.4 

306.2 

Lbm 

750 

4 18. 5 . 

422.3 

j 

429.8 

436.4 

442.3 

448.5 

455.6 

Ft 2 . See 

1000 

557.6 

561.2 

570.5 

579.0 

586.6 

594.7 

603.9 

Frozen 

100 

1.31 

1.31 

1.32 1 

1.33 

1.34 

1.34 

1.35 

Specific 

200 

1.31 

1.31 

1.32 

1.33 

1.33 

1.34 

1.34 

Heat 

300 

1.31 

1.31 

1.31 

1.32 

1.33 

1.33 

1.34 

Ratio 

500 

1.31 

1.30 

1.31 

1.32 

1.33 

1.33 

1.34 

7 FR 

750 

1.31 

1.30 

1.31 

1.32 

1.32 

1.33 

1.33 


1000 



1.30 

1.31 

1.31 

1.32 

1.32 

1.33 


NOTE: P c (psia) x 6.895 x 10 3 = P c (w/m 2 ) 


Mass Flow ( — — ) 
Ft .sec 


X 4.883 = MaSS FlOW ( ~ry ^ ) 

sec 




FIG. 1-1 
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VARIATION OF EQUILIBRIUM VACUUM SPECIFIC 
IMPULSE WITH AREA RATIO 

H 2 (1) - F 2 (i) 

P c = 100 PSIA (6.895X10 5 N/m2) 
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VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 


FIG. 1-2 
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VARIATION OF EQUILIBRIUM VACUUM SPECIFIC 
IMPULSE WITH AREA RATIO 


FIG. 1-5 


H 2 (fl - F 2 (i) 

P c = 300 PSIA |(2.069 X 10 6 N/m 2 ) 



EXIT AREA RATIO, (A/A M|N ) EX1T 





























H910666-16 FIG. 1-6 

VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE ilTH AREA RATIO 

H 2 (i) - P 2 CO 

P c =300 PSI A f(2,069 X 10 6 N/m 2 ) 
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VARIATION OF EQUILIBRIUM VACUUM SPECIFIC 
IMPULSE WITH AREA RATIO 


h 2 (i) - f 2 (l) 

P c = 500, P SI A | (3.448 X 10 6 N/m 2 ) 































FIG. 1-8 
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VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

H 2 CO - F 2 CO 

P c = 500 PSIA (3.448 X 10 6 N/m 2 ) 
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FIG. 1-9 


VARIATION OF EQUILIBRIUM VACUUM SPECIFIC 
IMPULSE WITH AREA RATIO 


H 2 W - F 2 (i) 

P c = 750 PSIA (5.171 X 10 6 N/m 2 ) 



EXIT AREA RATIO, (A/A M|N ) EX1T 
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VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

H 2 (0 - f 2 CO 

P C = 750 PSIA (5.171 X 10 6 N/m 2 ) 


FIG. 1-10 
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VARIATION OF EQUILIBRIUM VACUUM SPECIFIC 
IMPULSE WITH AREA RATIO 


FIG. 1-11 


H 2 (£> - F 2 (i) 

P c =11000 PSIA (6.895X10 6 N/m 2 ) 



EXIT AREA RATIO, (A/A M)N ) EX1T 


















H910666-16 FIG. 1-12 

VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

h 2 (j0 - f 2 CjO 

P c =1000PSIA (6.895X10 6 N/m 2 ) 
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FIG. 1-13 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


H 2 it) - F 2 (i) 

P c = 100 PSIA (6.895X105N/m 2) 

• 0/F = 8.0 

NOTE: REACTION RATE CONSTANTS LISTED IN TABLE 1-1 
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H9 10666— 16 FIG. 1-14 

NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

H 2 (f ) - F j (i) 

P c = 100 PSIA (6.895X10® N/m 2) 

0/F = 12,0 

NOTE: REACTION RATE CONSTANTS LISTED IN TABLE 1-1 








H ATOM GRADIENT, (8jH|/at) p - kg-MOLE m" 3 SEC" 
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FIG. 1-15 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

H 2 (!) - F j (i) 

P c = 100 PSIA(6.895X10 s N/m2) 

0/F = 16.0 

NOTE: REACTION RATE CONSTANTS LISTED IN TABLE I- 1 
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FIG. 1-16 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

H 2 <£) - F 2 (l) 

P c = 200PSIA(1.379X10<>N/m2). 

0/F = 8.0 

NOTE: REACTION RATE CONSTANTS LISTED IN TABLE 1-1 







H ATOM GRADIENT, (8[H|/8t)p - kg-MOLE m" 3 SEC" 
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FIG. 1-17 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 
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P c =200 PSI A ( 1.379X10 ^ N/m2) 

0/F = 12.0 

NOTE: REACTION RATE CONSTANTS LISTED IN TABLE. 1-1 
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FfG. 1-18 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


h 2 u ) - f 2 at) 

P c = 200 PSIA(1.379X106N/m2) 


0/F = 16.0 

NOTE: REACTION RATE CONSTANTS LISTED IN TABLE 1-1 
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.01 
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FIG. 1-19 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


H 2 (i ) - F j (£) 

P c = 300 p SI A (2.069X10 6 N/m 2 ) 
0/F = 8.0 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE 1-1 










H ATOM GRADIENT, (8[H]/8t) D - kg-MOLE m" 3 SEC 
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FIG. 1-20 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

H 2 (i ) - F 2 (i) 

P c = 300 PSIA (2.069X10 6 N/m 2 ) 

0/F = 12.0 

NOTE: REACTION RATE CONSTANTS LISTED IN TABLE 1-1 
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FIG. 1-21 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


H 2 U ) - F 2 
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P c = 300 PSIA (2.069X10 6 N/m 2 ) 

O/F = 16.0 

NOTE: REACTION RATE CONSTANTS LISTED IN TABLE 1-1 
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FIG. 1-22 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

h 2 a ) - f 2 a) 

P c = 500 PSIA (3.448X10 6 N/m 2 ) 

0/F = 8.0 

NOTE: REACTION RATE CONSTANTS LISTED IN TABLE 1-1 









H ATOM GRADIENT, OM/9t) p - kg-MOLE m~ 3 SEC" 
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NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 1-23 


H 2 U)- F j I*) 

P c = 500 PSIA (3.448X10 6 N/m 2 ) 

0/F = 12.0 

NOTE: REACTION RATE CONSTANTS LISTED IN TABLE 1-1 




H ATOM GRADIENT, (9[H|/9t) p - kg-MOLE m“ 3 SEC" 


1910666- 16 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 1-24 


h 2 <i ) - f 2 a) 

P c = 500 PSIA (3.448X10 6 N/m 2 ) 


0/F = 16.0 

NOTE: REACTION RATE CONSTANTS LISTED IN TABLE 1-1 
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NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 1-25 


H 2 <£ ) - F 2 (i) 

P c = 750 PSIA (5.171X10 6 N/m 2 ) 
0/F = 8.0 



SUPERSONIC AREA RATIO, A/Amin 


SUBSONIC 




H ATOM GRADIENT, 0|HJ/3t) D - kg-MOLE SEC 
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FIG. 1-26 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

H 2 (!) - F 2 (i) 

P c = 750 PSIA(5.171X10 6 N/m2) 

0/F = 12.0 

NOTE: REACTION RATE CONSTANTS LISTED IN TABLE 1-1 



SUBSONIC - U SUPERSONIC AREA RATIO, A/Amin 
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10666- 16 


FIG. 1-27 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

H 2 a ) - F j (I) 

P c = 750 PSIA (5.171X10 6 N/m 2) 


0/F = 1 &0 

NOTE: REACTION RATE CONSTANTS LISTED IN TABLE 1-1 



SUBSONIC — *-l»* — SUPERSONIC AREA RATIO, A/Amin 




H ATOM GRADIENT, 0[H]/9t) p - kg-MOLE m" 3 SEC" 
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FIG. 1-28 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


H 2 (i ) - F 2 U) 

P c =1000 PSIA (6.895X10 6 N/m 2 ) 
0/F = 8.0 


10,000 


1000 


100 


1.0 





H ATOM GRADIENT, OlH]/at) p - kg-MOLE m “ 3 SEC" 
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FIG. 1-29 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


h 2 m - f 2 (i) 

P c =1000 PS1A (6.895X106N/m2) 


0/F = 1Z0 

:: REACTION RATE CONSTANTS LISTED IN TABLE 1-1 


COMPOSITE KINETIC I 
(TRANSITION FACTOR = 2 - 


TABLE OF REACTIONS 


REACTION 


HfH+F - 
H + F + F - 
H + H + H - 

H + F + H ■ 
H + H+HF* 
H + F -f HF * 

H + H + H 2 - 
H + F + Ho - 


H 2 + F 
HF + F 
H 2 fH 
HFtH 
H 2 + HF 
HF + HF 
H 2 + H 2 
HF + H 2 


EQUILIBRIUM ( 15 DEG CONE, r f = 1.0 FT(0. 3048m) 

ITT 1 ITTl T f c /r t = H+HT+-H-H-T H ~f 

\======4# i # h ======~=“ i; 

_a_ y=|EQUI LIBRIUM || II I II I ' 

~~~\T -0/r, ) [d(A / A M|N ) / d(x/r,)]=lF 

J. . = 1#0 FT (0.3048 m)HTTH+HT 


IlSIJISIIISi 


I 



saSS 
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■ii 



»<■■■■■■■■■■■■! 

liSIUIIIHU! 
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H 



mm 


SUBSONIC- 


■ SUPERSONIC AREA RATIO, A/Amin 



H ATOM GRADIENT, OlH]/9t) p - kg-MOLE m~ 3 SEC - 
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FIG. 1-30 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 
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’ c =1000 PSIA (6.895X10 6 N/m 2 ) 
0/F = 16.0 
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FIG. 1-31 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-FLUORINE PROPELLANT SYSTEM 


H 2 (i) - F 2 (i) 

P c = 100 PS1A (6.895X10 5 N/m2) 
0/F = 8.0 



FREEZING AREA RATIO, (A/A M1N ) FR 
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FIG. 1-32 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN PROPELLANT SYSTEM 


H 2 (!) - F 2 (£) 

P c = 100 P SI (6.895X10 5 N/m2) 
0/F = 12.0 



FREEZING AREA RATIO, (A/A MIH ) FR 
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FIG. 1-33 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-FLUORINE PROPELLANT SYSTEM 
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FIG. 1-34 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-FLUORINE PROPELLANT SYSTEM 


H 2 (i) - F 2 U) 

R C = 200 PSIA (1.379X10 6 N/m 2 ) 
O/F = 8.0 
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FIG. 1-35 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-FLUORINE PROPELLANT SYSTEM 

Hj (!) - F 2 (!) 

P c = 200 PSIA (1.379X10 6 N/m 2 3 ) 

0/F = 12.0 


O 


iu 


u 

< 

> 

cl 


•jt 


m 

%A 

-j 

a. 

2 

y 

C 

u 

UJ 

a. 

to 


2 

3 

P 

U 

< 

> 


4800 


4700 


4600 


4500 


4400 


4300 


4200 


4100 
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FIG. 1-36 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-FLUORINE PROPELLANT SYSTEM 


H 2 (4) - f 2 U) 

P c = 200 PSIA (1.379X10 6 N/m 2 ) 
0/F = 16.0 



FREEZING AREA RATIO, (A/A MIN ) FR 
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FIG, 1-37 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-FLUORINE PROPELLANT SYSTEM 


H 2 (£) - F 2 (£) 

P C = 300 PSIA (2.069 X 10 6 N/m 2 ) 
0/F = 8.0 



FREEZING AREA RATIO, (A/A MIN ) FR 
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FIG. 1-38 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-FLUORINE PROPELLANT SYSTEM 

H 2 (£) - F 2 (l) 

P c = 300 P SI A (Z069 X 10 6 N/m 2 ) 

0/F = 12.0 
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FIG. 1-39 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-FLUORINE PROPELLANT SYSTEM 


H 2 ($) - F 2 (?) 

P c = 300 PSIA (2.069 X 10 6 N/m 2 ) 
0/F = 16.0 
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FIG. 1-40 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-FLUORINE PROPELLANT SYSTEM 

H 2 Cf) - F 2 (jD 

P c = 500 PSIA (3.448 X ]0 6 N/m 2 ) 

0/F = 8.0 



FREEZING AREA RATIO, (A/A M1N ) FR 
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FIG. 1-41 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-FLUORINE PROPELLANT SYSTEM 


h 2 it) - f 2 (d 

P c = 500 PSIA (3.448 X 10 6 N/m 2 ) 
0/F = 12.0 



FREEZING AREA RATIO, (A/A MIN ) FR 
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FIG. 1-42 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-FLUORINE PROPELLANT SYSTEM 


h 2 a) - f 2 to 

P c = 500 PSIA (3.448 X 10 6 N/m 2 ) 
0/F = 16.0 



FREEZING AREA RATIO, (A/A min ) F r 
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FIG. 1-43 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-FLUORINE PROPELLANT SYSTEM 

H a (0 - F 2 (I) 

P c = 750 PSIA (5.171 X 10 6 N/m 2 ) 

0/F = 8.0 
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FIG. 1-44 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-FLUORINE PROPELLANT SYSTEM 
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P c = 7 50 PS IA (5.171 X 10 6 N/m 2 ) 
O/F = 12.0 
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FIG. 1-45 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-FLUORINE PROPELLANT SYSTEM 


h 2 U) - F 2 (0 

P c = 750 P5IA (5.171 x 10 6 N/m 2 ) 


0/F = 16.0 



FREEZING AREA RATIO, (A/A MIN ) FR 
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FIG. 1-46 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-FLUORINE PROPELLANT SYSTEM 


H 2 (I) - F 2 U) 

P c =1000PSIA (6.895X106N/m2) 
0/F = 84) 



FREEZING AREA RATIO, (A/A MIN ) FR 
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FIG. 1-47 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-FLUORINE PROPELLANT SYSTEM 


h 2 a) - f 2 to 

P c =1000PSIA (6.895X10 6 N/m 2 ) 
0/F = 12.0 



FREEZING AREA RATIO, (A/A MIN ) FR 
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FIG. 1-48 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-FLUORINE PROPELLANT SYSTEM 


H 2 (i) - F 2 (/) 

P c =1000 P SI A (6.895X10 6 N/m 2 ) 
0/F =16.0 
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TABLE II -1 

SUMMARY OF ELEMENTARY REACTIONS AND REACTION RATE 
CONSTANTS EMPLOYED IN H 2 -0 2 RECOMBINATION MECHANSIM 


REACTIONS 

FORWARD RATE 


H + H + A — 

H 2 + A 

kf = 4.62 X 10 14 T“ 1 


H + OH + A — 

H 2 0 + A 

k { = 7.85 X 10 15 T“\ 

1 

H+H + H 

H H 

k f = 20 (4.62 X 10 14 T“ 1 ) 

2 

h + h + h 2 — 

H . + H „ 
2 2 

k f = 2.5 (4.62 X lO^T- 1 ) 

3 

h + h + h 2 o 

h 2 + h 2 o 

K f = 2.5 (4.62 X 10 15 T -1 ) 

4 

h + oh + h 2 o — 

h 2 o + H 2 0 

k, = 20 (7.85 X 10 15 T -1 ) 

5 

H 4- OH + Ho — ^ 

h 2 o +-h 2 

k f = 3 (7.85 X 10 15 T* 1 ) 

6 

H + OH + H 

h 2 o + h 

k { = 3 (7.85 X I0 15 T -1 } 


NOTE: ALL RATES ARE EXPRESSED IN TERMS OF LB-MOLES, FT 3 , SEC; AND T IS IN °R 


TO CONVERT TO kg-MOLES, m 3 , SEC, AND T IN ° K; MULTIPLY ABOVE RATE BY 2.165 X UT 3 
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TABLE II-2 


EQUILIBRIUM AND FROZEN FLOW PROPERTIES AT 
NOZZLE THROAT' 

*2 CO - 0 2 (i) 



p c -psia 

OXIDIZER - FUEL WEIGHT RATIO 

3.0 

4.0 

5.0 

6.0 

7.0 

8.0 

Equilibrium 
Mass 
Flow - 

w/a 

Lba 

Ft 2 . Sec 

100 

58.2 

58.7 

60.4 

62.6 

65.0 

67.5 

200 

116.5 

117.1 

120.1 

124.2 

129.0 

133.9 

300 

174.7 

175.5 

179.6 

185.6 

192.6 

199.9 

500 

291.1 

292.2 

298.4 

307.8 

319.2 

331.2 

750 

436.7 

437.9 

446.6 

460.1 

476.6 

494.6 

1000 

582.2 

583.7 

594.6 

611.9 

633.6 

657.4 








Frozen 

Mass 

Flow - 

W/A 

Lbm 

Ft 2 . Sec 

1 

100 

58.6 

59.7 

61.8 

64.3 

66.9 

69.4 

200 

117.0 

118.8 

122.6 

•127.3 

132.4 

137.4 

300 

175.3 

177.7 

183.1 

190.0 

197.4 

204. 9 

500 

292.0 

295.2 

303.6 

314.4 

326.5 

339.0 

750 

437.7 

441.9 

453.6 

469.3 

487.0 

505.5 

1000 


588.5 

603.3 

623.5 

646.8 

671.3 








Frozen 

Specific 

Heat 

Ratio 

Vfr 

100 

1.25 

1.23 

1.22 

1.21 

1.21 

1.21 

200 

1.25 

1.23 

1.22 

1.21 

1.21 

1.21 

300 

1.25 

1.23 

1.22 

1.21 

1.20 

1.20 

500 

1.25 

1.23 

1.21 

1.21 

1.20 

1.20 

750 

1.25 

1.23 

1.21 

1.20 

1.20 

1.20 

— 

1.25 

1.23 

1.21 

1.20 

1.20 

1.20 


NOTE: P c (psia) x 6.^95 x K>3 = p c {N/m 2 ) 

Mass Flow ( ^bm ) x 4.883 = Mass Flow ( , T k § — ) 
Ft2. sec 

























































































































H910666- 16 FIG. 11-1 

VARIATION OF EQUILIBRIUM VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

h 2 <i)'-o 2 U) 

P c = 100 PSIA (6.895X105N/m2) 







































H9 10666- 16 FIG. 11-2 

VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 


H 2 (i) -0 2 (i) 

P c = 100 PSIA (6.8WX10 5 N/m2) 
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H9 10666- 16 

VARIATION OF FROZEN VACUUM. SPECIFIC IMPULSE WITH AREA RATIO 

H 2 (!) -0 2 (!) 

P c = 200 PSIA (1.379X106N/m2) 






















































FIG. 11-6 
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VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

H 2 (J) - 0 2 (|) 

P C =300 PSIA (2.069X106N/m2) 


4500 


4400 


O) 

S 4300 
u 

iu 

«/> 


z 

I 

u 4200 

> 

0. 


IU 

•J 

a. 

35 


± 4100 


u 

in 

a. 

07 

X 

3 

=> 

U 

< 

> 


4000 


3900 


3800 













H9 10666- 16 FIG. 11-7 

VARIATION OF EQUILIBRIUM VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

h 2 <i) - o 2 (l) 

P c = 500 PSIA (3.448X106N/ m 2) 
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H9 10666- 16 FIG. 11-8 

1 VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

h 2 <i) - o 2 (i) 

P c = 500 PSIA ,(3.448X10 6 N/m2), 
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VACUUM SPECIFIC IMPULSE, T SP - (N)(SEC)/(kg) 

SP VAC 
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VARIATION OF EQUILIBRIUM VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

H 2 (i) - 0 2 (M 

P c = 7 50 PSIA (5.171X10 6 N/m 2 ) 
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VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

h 2 U) -o 2 (i) 

P c = 750 PSIA (5.171X106N/m2) 


4500 


4400 


u 

IU 

CO 


4300 


z 

i 


% 4200 
> 




Ul 

V) 


CL 

2 4100 

u 


u 

UJ 

CL 

%/> 


U 

< 

> 


4000 


3800 


3800 






VACUUM SPECIFIC IMPULSE, T^ p - (N)(SEC)/(kg) 


H9 10666- 16 FIG. 11-11 - 

VARIATION OF EQUILIBRIUM VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

H 2 (£) - 0 2 (£) 

P c =1000 PSIA (6.895X106N/m 2) 



0 20 40 60 80 100 

EXIT AREA RATIO, (A/A MIN ) 











VACUUM SPECIFIC IMPULSE, I, p - (N) (SEC)/(kg) 

5t VAC 


H9 10666- 16 FIG. 11-12 

VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

H 2 (£) - 0 2 (£) 


P C =1000PSIA ,(6.895X10 6 N/m 3 



0 20 40 60 80 100 

EXIT AREA RATIO, <A/A M|N ) Em 















H ATOM GRADIENT, 0[H] /8*)_ - kg-MOLE nr 3 SEC" 


H9 10666- 16 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 11-13 


H j (i ) — 0 j (i) 

P c = 100 PSIA (6.895X105N/m2) 
0/F = 4.00 



2 1 2 3 4 ' 5 6 

SUBSONIC' » - l— SUPERSONIC AREA RATIO, A/A M|N 




H ATOM GRADIENT, 0[H] /3t)„ - kg-MOLE m~ 3 SEC~ 
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NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 11-14 


H 2 U)-0 2 (l) 

P c = 100 PSIA (6.895X105N/m2) 
0/F = 6.00 
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NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 11-15 


H j (i) - O j (£) 

P c = 100 PSIA (6.895 X 10 5 N/m 2 ) 
O/F = 8.00 



subsonic— SUPERSONIC AREA RATIO, A/A M1N 


H ATOM GRADIENT, ON) /d*) 0 - kg-MOLE m" 3 SEC” 


H9 10666- 16 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 11-16 


H j (i ) - 0 2 U) 

P c = 200 PSIA ;(1.379X106N/m 2) 
O/F = 4.00 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE !||-1 





H ATOM GRADIENT, 0[H] / 9 ») - kg-MOLE m" 3 SEC" 
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NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 11-17 


H j (i) - 0 j (1) 

P c = 200 PSIA (1.379X106N/m 2) 
0/F = 6.00 
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FIG. 11-18 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

H j (i) - 0 j (i) 

P c = 200 PSIA (1.379X106N/m2) 

0/F = 8.00 



H ATOM GRADIENT, 0|H] /9tL - kg -MOLE m~ 3 SEC" 
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NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 11-19 


h 2 (|)-o 2 (1) 

P c = 300 PSIA (2.069X10 6 N/m 2 ) 
0/F = 4.00 



SUBSONIC-—!— —SUPERSONIC AREA RATIO, A/A M)N 
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NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 11-20 


H 2 (i)-0 2 (JE) 

P c = 300 PSIA (2.069X10 6 N/tn 2 ) 
O/F = 6.00 



SUBSONIC- -*^-* — SUPERSONIC AREA RATIO, A/ A M|N 
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NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 11-21 


H 2 (£) - 0 2 U) 

P c = 300 PSIA(2.069X10 6 N/m2) 
0/F = 8.00 
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FIG. 11-22 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

H j (i ) - 0 j (1) 

P c = 500 PSIA (3.448X10 6 N/m2) 

0/F = 4.00 



SUBSONIC— *-}-«— : SUPERSONIC AREA RATIO, A/A M(N 




H ATOM GRADIENT, 0[H] /9t)_ - kg-MOLE m 





SUBSONIC 


SUPERSONIC AREA RATIO, A/ A. 







H ATOM GRADIENT, 0|H] /dt) - kg-MOLE m~ 3 SEC" 


FIG. 11-24 
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NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


H 2 (i) - 0 2 (£> 

P c = 500 PSIA (3.4'48X10 6 N/m 2 ) 
0/F = 8.00 




H ATOM GRADIENT, 0|H] /9») n - kg-MOLE m~ 3 SEC“ 
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NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 11-25 


H j (£) - 0 2 (i) 

P c = 750 P SI A (5.171X10 6 N/m2) 
0/F = 4.00 
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NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 11-26 


H 2 (i)-0 2 (£) 

P c = 750 PSIA (5.171X106N/m2) 
0/F = 6.00 



SUBSONIC — SUPERSONIC AREA RATIO, A/A M|N 






H ATOM GRADIENT, 0[H)/9») p - kg-MQLE m~ 3 SEC” 


H9 10666- 16 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 11-27 


H 2 U)- 0 2 U) 

P c = 750 PSIA (5.171X106N/m2) 
0/F = 8.00 



NOTE; REACTION RATE CONSTANTS LISTED IN TABLE II — 1 



? . 1 2 3 4 5 6 


SUBSONIC 
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AREA RATIO, A/A M|N 





H ATOM GRADIENT, (9[H|/9t) p - kg-MOLE m“ 3 SEC” 
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NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 11-28 


H j (i ) - O j (i) 

P c =1000 P SI A (6.895X106 N/m2) 
0/F = 4.00 




H ATOM GRADIENT, (8[H|/8t) p - kg-MOLE m“ 3 SEC” 


H9 10666- 16 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 11-29 


H 2 (i) - 0 j U) 

P c =1000 PSIA (6.895X10 6 N/m2) 
0/F = 6.00 



2-1 2 3 4 5 6 

SUBSONIC - * ^— 1 - SUPERSONIC AREA RATIO, A/A M|N 




H ATOM GRADIENT, 0[H]/8t) p - kg-MOLE SEC” 
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NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 11-30 


H j (&) “ 0 j (i) 

P c =1000 P SI A (6.895X10 6 N/m 2 ) 
0/F = 8.00 



■j** — SUPERSONIC AREA RATIO, A/ A m|n 


SUBSONIC 
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FIG. 11-31 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-OXYGEN PROPELLANT SYSTEM 

H j (i) -0 2 U) 

P c = 100 PSIA (6.895X105N/m2) 

0/F = 4.00 



FREEZING AREA RATIO, <A/A M1N ) pR 
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FIG. 11-33 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-OXYGEN PROPELLANT SYSTEM 

H 2 <!) - Q 2 U) 

P c = 100 PS1A (6.895X105N/m 2) 

0/F = 8.00 
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FIG. 11-36 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-OXYGEN PROPELLANT SYSTEM 

H 2 (f) -0 2 (i) 

P c = 200 P$IA,(1.379X106N/m2) 

0/F = 8.00 
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FIG, 11-37 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-OXYGEN PROPELLANT SYSTEM 

H 2 (i) -0 2 (i) 

P c = 300 PSIA ( 2 . 069 X 106 N/m 2 ) 

0 /F = 4.00 
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FIG. 11-38 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-OXYGEN PROPELLANT SYSTEM 

H 2 (i) -0 2 (jH 

P c = 300 PSIA i(2,069X10 6 N/m2) 

0/F = 6.00 
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FREEZING AREA RATIO, (A/A M1N ) pR 
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FIG. 11-40 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-OXYGEN PROPELLANT SYSTEM 

H 2 (i) - 0 2 (i) 

P c = 500 PSIAi(3.448X10 6 N/m 2 ) 

0/F = 4.00 
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FIG. 11-41 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-OXYGEN PROPELLANT SYSTEM 


H j (i ) -0 2 (i) 

P c = 500 PSIA (3.448X106N/m2) 
0/F = 6.00 
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FIG. 11-42 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-OXYGEN PROPELLANT SYSTEM 

H 2 (f ) - 0 2 (i) 

P c = 500 PSIA (3.448X106N/m 2 ) 

0/F = 8.00 
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FIG. 11-44 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-OXYGEN PROPELLANT SYSTEM 


H 2 (iJ -0 2 (jO 

P c = 750 PSIA (5. 171X10 6 N/m2) 
0/F = 6.00 
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FIG. 11-45 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-OXYGEN PROPELLANT SYSTEM 


H j (i ) -0 2 Ut) 

P c = 750 PSIA(5.171X10 6 N/m2) 
0/F= 8.00 
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FJG. 11-46 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-OXYGEN PROPELLANT SYSTEM 


H 2 (i) -0 2 U) 

P c =1000 PSIA (6.895XlQ 6 N/m 2) 
0/F = 4.00 
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FIG. 11-47 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-OXYGEN PROPELLANT SYSTEM 


H 2 (|) -0 2 (i) 

P c =1000PSIA(6.895X106N/m2) 

0/F = 6.00 


< A/ ' A MIN * 
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FIG. 11-48 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR HYDROGEN-OXYGEN PROPELLANT SYSTEM 


H 2 m -o 2 u) 

P c =1000 PSIA (6.895X106 N/m 2) 
0/F = 8.00 
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TABLE III -1 

SUMMARY OF ELEMENTARY REACTIONS AND REACTION RATE CONSTANTS 
EMPLOYED IN CH 4 - FLOX RECOMBINATION MECHANISM 


REACTION 

FORWARD RATE 


H + H + H 


— 

H 2+H 

k f =4.62 X 10“ 54 T 

1 

H + H + H 


H 2 h- H 

k f = 20 (4.62 X 10 -14 T -1 ) 

2 



m mm 

kf = 2.5 (4.62 X 10~ 14 T -1 ) 

3 

H +H + HF 

~^z 

H 2 + HF 

k f = 2.5 (4.62 X 10-14 T' 1 ) 

4 

H + H + C0 2 


H 2+ co 2 

k f = 2.0(4.62 X 10“ 14 T _1 ) 

5 

H + H + CO 


h 2 + CO 

k f = ZO (4.62 X 10 -14 T~ 1 ) 

6 

H + H + H 2 0 


h 2 + h 2 o 

k f = 2.5 (4.62 X 10~ 14 T~’) 


H + F + M 


HF ^ M 

kf = 1.155 X 10 15 T _1 

7 

H + F + H 

— 

HF + H 

kf = 1 (1.155 X 10 15 T" 1 ) 

8 

H + F + H 2 


HF + H 2 

k, =2.5 (1.155 X 10 1 5 T -1 ) 

9 

H + F + HF 


HF + HF 

kf = 2.5(1.155 X 10 15 T _1 ) 

10 

H + F + C0 2 


HF + C0 2 

kf = ZO (1.155 X 10 15 T ~S 

n 

H + F + CO 


HF + CO 

kf = 2.0 (1.155 X 10 15 T _1 ) 

12 

H + F + H 2 0 


h f + h 2 o 

kf = 2.5(1.155 X 10 15 T -1 ) 


H + OH + M 


H 2 0 +M 

kf =7.85 X lO 1 ^ T _1 

13 

H + OH + H 


h 2 o + h 

kf = 3 (7.85 X 10 15 T -1 ) 

14 

h + oh + h 2 


h 2 o + h 2 

k f = 3 (7.85 X 10 15 T _1 ) 

15 

H + OH + HF 


h 2 o + hf 

kf =3 (7.85 X 10 1S T - 1 ) 

16 

H + OH + C0 2 


h 2 o + co 2 

kf = 10 (7.85 X 10 15 T- 1 ) 

17 

H + OH + CO 


h 2 o + CO 

k, =3(7.85 X 10 ,S T~’ ) 

18 

H +OH +H 2 0 


H 2 0 + H 2 0 

k, = 20 (7.85 X 10 15 T ~ 1 ) 


NOTE: ALL RATES ARE EXPRESSED IN TERMS OF LB-MOLES, FT 3 , SEC; AND T IS IN ° R 


TO CONVERT TO kg-MOLES, m? , SEC, AND T IN ° K; MULTIPLY ABOVE RATE BY 2.165 X 10~ 3 
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TABLE III -2 

EQUILIBRIUM AID FROZEN FLOW PROPERTIES AT 
NOZZLE THROAT 
CH^ - FLOX (82.6/0 Fg) 



P G -psia 

OXIDIZER - FUEL WEIGHT RATIO 

3-5 

4.0 

4.5 

5.0 

5.5 

6.0 

6.5 

Equilibrium 

Mass 

Flow - 
W/A 

Lbm 

Ft 2 . Sec 

100 

mm 

70.0 

m 

1 

I 

67.1 

II 

200 

mm 

■91' 

1 

1 

132.0 


134.6 

300 

205.5 

r~ 

207.4 

M 

197.9 

■91 

198.0 

200.8 

500 

H88H 

■m 

E HI 

191 

■91 

E9H 

m 

750 

512.9 

513.9 

516.5 

489.1 

486.2 

488.8 

496.0 

1000 

682.5 

667.7 


Hi 

645.8 

649.2 

mm 




Hi 

H 

H 

H 

■Hi 

Frozen 

Mass 

Flow - 

W/A 

Lbrn 

100 

■HI 

71.2 

II 

■EH 

WRSEk 

n 

■99 

200 

142.8 

l 4 o .9 

139.6 

138.9 

138.6 

— 

l 4 l .4 

300 

213.0 

210.0 

208.0 

206.8 

206.4 

207.8 

210.6 

500 

352.6 

347.3 

343.8 

341.6 

340.8 

343.0 . 

347.8 

KJ ILL 

Ft 2 . Sec 

750 

526.3 

517.9 

512.4 

508.9 

507.5 

510.7 

■1111 

1000 

699.3 

687.8 

680.2 

675.3 

673.2 

677.4 - 

687.1 









Frozen 
Specific 
Heat 
Ratio 
V FR 

100 

■19 

1.30 

■EE 

■ m 

mm 

IRS 

1.34 

200 

1.28 

1.30 

1.31 

1.32 

1.32 

— 

■H 

300 

1.28 

1.29 

1.30 

1.31 

1.32 

H&IM 

HH 

mm 

1.28 

mm 

1.30 

1.31 

1.32 

1.32 _ 

II 

■n 

1.28 

1.29 

1.30 

mm 

1.32 

_ 1.32 

1.32 

1000 

1.28 

1.29 

1.30 

1.30 

mmmm 

Wmmm 

1.32 

1.32 


NOTE: P G (psia) x 6.895 x 103 = p (N/m 2 ) 

Mass Flow ( — ) x 4.883 = Mass Flow ( -n - ~ ) 

Ft^.sec rsr sec 
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FIG. 111-1 


VARIATION OF EQUILIBRIUM VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

CH 4 - FLOX (82.6% F 2 ) 

P c = 100 PS1A (6.895 X 10 5 N/m 2 ) 
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FIG. 111-2 


VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 


CH 4 - FLOX (82.6% F 2 ) 

P c = 100 PSIA (6.895 X 10 5 N/m 2 ) 



0 20 40 60 80 100 

EXIT AREA RATIO, (A/A M1N ) EXjT 
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FIG. 111-3 


VARIATION OF EQUILIBRIUM VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

CH 4 - FLOX (82.6% ? 2 ) 

P c = 200 PSIA (1.379 X 10 6 N/m 2 ) 
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FIG. 111—4 


VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 


CH 4 - FLOX (82.6% ) 

P c = 200 PSI A ( 1.379 X 10 6 N/m 2 ) 
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FIG. 111-5 


VARIATION OF EQUILIBRIUM VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

CH 4 - FLOX (82.6% F 2 ) 

P c =300 PSIA (2.069 X 10 6 N/m 2 ) 
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FIG. 111-6 


VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

CH 4 - FLOX (82.6% F 2 ) 

P c = 300 PSIA (2.069 X 10 6 N/m 2 ) 
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FIG. Ill— 7 


VARIATION OF EQUILIBRIUM VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

CH 4 - FLOX (82.6% F 2 ) 

P c =500 PSIA (3.448 X 10 6 N/m 2 ) 
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FIG. I II -8 


VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

CH 4 - FLOX (82.6% F 2 ) 

P c =500 PSIA (3.448 X 10 6 N/m 2 ) 
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FIG. 111-9 


VARIATION OF EQUILIBRIUM VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

CH 4 - FLOX (82.6% ) 

P c =750 PSIA (5.171 X 10 6 N/m 2 ) 
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FIG. 111-10 


VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 


CH 4 - FLOX (82.6% F 2 ) 

P c = 750 PSIA (5.171 X 10 6 N/m 2 ) 
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FIG. 111-11 
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FIG. 111-12 


VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

CH 4 - FLOX (82.6% F,, ) 

P c =1000 PS1 A (6.895 X 10 6 N/m 2 ) 
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FIG. 111-13 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


CH 4 - FLOX (82.6% F 2 ) 

P c = 100 PSIA (6.895X10 5 N/m2) 



2 1 2 3 4 5 6 

f— SUPERSONIC AREA RATIO, A/ A m , n 


SUBSONIC 




H ATOM GRADIENT, OM/9») p - kg-MOLE m~ 3 SEC” 
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FIG. 111-14 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


REACTION 

19, 

1 ST 


CH 4 - FLOX (82.6% F 2 ) 

P c = 100 PSIA (6.895X10^N/m 2) 

O/F - 5.50 

RATE CONSTANTS LISTED IN TABLE |||-1 


TABLE OF REACTIONS 


-EQUI LIBRIUM ~j— p 1 1 1 1 H | M 0 
-0/r t ) (A/A m|n ) / d(x/f f )]= 7 

" 1 FT" 1 r = i.o FT (0.3048 m) ~j ® 


REACTION 


H + H + H 
H + H + H 2 
H + H + HF 
H 4 H + C0 2 
H + H + CO 
H + H . + H 2 0 
H + F + H 

h + f + h 2 

H+F+HF 
H + F + C0 2 
H + F + CO 
H + F + H 2 0 

H+OH+H 
H + OH + H 2 
H+OH+HF 


H 2 + H 

■h 2 + h 2 

H 2 +HF 

-h 2 + co 2 

•h 2 + CO 
h 2 + h 2 o 

HF +H 
HF + H 2 
HF + HF 
HF + C0 2 

HF + CO 
-HF + H 2 0 

- H 2 0 + H 
■ H 2 0 + H 2 

-H 2 0 + HF 


H + OH + C0 2 ->— H 2 0 + C0 2 


H + OH + CO 


-h 2 o + CO 


H + OH + H 2 O^H 2 0 + H 2 0 


♦EQUILIBRIUM (15 DEG CONE, r f = - 

J FT (0.3048 m)[, r~/7 
a 1 1 111 1 1 ji 1 i 1 1 1 r 1 1 

Sj-COMPOSITE KINETIC 

1 (TRANSITION FACTOR = 2) 



-COMPOSITE KINETIC 
(TRANSITION 
‘^□FACTOR = 1) 


mm, 


iisWmisi 



SUBSONIC- 


-SUPERSONIC AREA RATIO, A/ A M1N 


H ATOM GRADIENT, OM/8*) p - kg-MOLE m" 3 SEC~ 
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NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG, 111-15 


CH 4 - FLOX (82.6% F 2 ) 

P c = 100 PSIA (6.895X10 5 N/m 2) 
0/F = 6.50 
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FIG. Ill -16 

NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


CH 4 - FLOX (82.6% F 2 ) 

P c = 200 PSIA (1.379X106N/m2) 
0/F = 3.50 



SUBS0NIC-*4^-SUPERS0NIC AREA RATIO, A/A min 


H ATOM GRADIENT, 0[Hj/at) p - kg-MOLE m~ 3 SEC* 


FIG. 111-17 
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NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


CH 4 - FLOX (82.6% F 2 ) 

P c = 200 PSIA (1.379X106N/m2) 
0/F = 5.50 



SUPERSONIC AREA RATIO, A/A M|N 


SUBSONIC 



H ATOM GRADIENT, 0[H|/8t) p - kg-MOLE m~ 3 SEC* 
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NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 111-18 


CH 4 - FLOX (82.6% F 2 ) 

P c = 200 PSIA (1.379X10 6 N/m 2) 
0/F = 6.50 



SUBSONIC— PE RSONIC AREA RATIO, A/A M)N 



H ATOM GRADIENT, OM/9t) p - kg-MOLE m' 3 SEC“ 
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<1 


FIG. 111-19 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


CH 4 - FLOX (82.6% F 2 ) 


P c = 300 PSIA (2.069X10 ^N/m 2) 



2 1 2 3 4 5 6 

SUPERSONIC AREA RATIO, A/A M1N 


SUBSONIC 



H ATOM GRADIENT, 0[Hj/9t) p - kj-MOLE m~ 3 SEC~ 
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FIG. 111-20 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

CH 4 - FLOX (82.6% F 2 ) 

P c = 300 PSIA (2.069X106 N/m 2) 

O/F = 5.50 





H ATOM GRADIENT, (8fH|/8t)p - k«-MOLE m" 3 SEC* 
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FIG. 111-21 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

CH 4 - FLOX (82.6% F 2 ) 

P c = 300 PSIA (2.069X10 6N/m2) 

0/F - 6.50 



SUBSONIC — SUPERSONIC AREA RATIO, A/A M|N 



H ATOM GRADIENT, 0[H|/9t) p - kg-MOLE m" 3 SEC" 
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FIG. 111-22 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


CH 4 - FLOX (82.6% F 2 ) 

P c = 500 PSIA(3.448X10 6 N/m2) 
0/F =3.50 



2 1 2 3 4 5 6 

subsonic-*4**-supersonic AREA RATIO, A/A M1N 


H ATOM GRADIENT, OlHj/8t) p - kg-MOLE m" 3 SEC* 
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FIG. 111-23 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


CH 4 - FLOX (82.6% F 2 ) 

P c = 500 PSIA (3.448X10^N/m 2) 
O/F - 5.50 



SUBSONIC-»>}-*-SUPERSONIC AREA RATIO, A/A m|N 


H ATOM GRADIENT, 0[H|/8t)p - kg-MOLE m~ 3 SEC" 
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FIG. 111-24 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


CH 4 - FLOX (82.6% ? 2 ) 

P c = 500 PSIA(3.448X106N/m2) 
0/F = 6.50 



SUBSONIC— SUPERSONIC AREA RATIO, A/A M|N 


H ATOM GRADIENT, 0[H}/9t) p - kg-MOLE m" 3 SEC~ 
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FIG. 111-25 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


CH 4 - FLOX (82.6% F 2 ) 

P c = 750 PSIA(5.171X106N/m2) 
O/F = 3.50 




H ATOM GRADIENT, 0[Hj/8t) p - kg-MOLE m" 3 SEC" 
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FIG. 111-26 

NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


CH 4 - FLOX (82.6% F 2 ) 

P c = 750 PSIA (5. 171X10^N/m 2) 
0/F = 5.50 
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NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 111-27 


CH 4 - FLOX (82.6% F 2 ) 
P c = 750 PSIA (5.171X10 6 N/m 2 ) 
0/F = 6.50 



SUBSONIC 




H ATOM GRADIENT, (3(Hj/9t) p - kg-MOLE m~ 3 SEC' 
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NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 111-28 


CH 4 - FLOX (82.6% F 2 ) 

P c = 1000 PSIA(6.895X106N/m 2 ) 
O/F = 4.00 










































H ATOM GRADIENT. (3|H|/at) p - k fl -MOLE m - 3 SEC“ 
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FIG. 111-29 

NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


CH 4 - FLOX (82.6% F 2 ) 
P c = 1000 PSIA (6.895X10 6 N/m2) 



SUBSONIC— *■]-*— SUPERSONIC AREA RATIO, A/A M|N 




H ATOM GRADIENT, 0(H|/at) p - kg-MOLE m~ 3 SEC " 
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NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


FIG. 111-30 


CH 4 - FLOX (82.6% F 2 ) 
P c = 1000 PS1A (6.895X106 N/m2) 
0/F - 6.50 



SUBSONIC — SUPERSONIC AREA RATIO, A/A M)N 
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FIG. 111-31 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR METHANE - FLOX PROPELLANT SYSTEM 

CH 4 - FLOX (82.6% ^ ) 

P c = 100 PSIA (6.895 X 10 5 N/m 2 ) 

0/F = 3.50 
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FIG. 111-32 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR METHANE - FLOX PROPELLANT SYSTEM 

CH - FLOX <82.6% F„ ) 

4 l 

P c = 100 PSIA <6.89 5 X 10 5 N/m 2 ) 

0/F=5.5 
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FIG. 111-33 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR METHANE - FLOX PROPELLANT SYSTEM 

CH 4 - FLOX (82.6% F, ■ ) 

P c = 100 PSIA (6.895X105 N/m2) 

0/F = 6.50 
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FIG. 111-34 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR METHANE - FLOX PROPELLANT SYSTEM 

CH 4 - FLOX (82.6% F 2 ) 

P c = 200 PSIA (1.379 X I0 6 N/m 2 ) 

0/F = 3.50 
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FIG. 111-35 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR METHANE - FLOX PROPELLANT SYSTEM 

CH 4 - FLOX (82.6% ^ ) 

P c = 200 PSIA ( 1.379 X 10 6 N/m 2 ) 

0/F = 5.50 


U 

UJ 

UO 

2 

I 


U 

< 

> 

H. 

i/> 

M 

uT 


3 

Ql 

U 

a. 



3 

U 

> 




H9 10666- 16 


FIG. 111-36 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR METHANE - FLOX PROPELLANT SYSTEM 

CH - FLOX (82.6% F ) 

4 2 

P C = 200 PSIA (1.379X106N/m2) 

0/F =6.50 
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FIG. 111-37 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR METHANE - FLOX PROPELLANT SYSTEM 

CH - FLOX '{82.6% F, ) 

4 2 

P c = 300 PSIA( 2.069 X 10 6 N/m 2 ) 

0/F = 3.50 
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FIG. 111-38 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR METHANE - FLOX PROPELLANT SYSTEM 


CH - FLOX (82,6% F ) 

4 2 

P c =300 PSIA (2.069 X 10 6 N/m 2 ) 
0/F = 5.50 



FREEZING AREA RATIO, (A/A MIN ) pR 
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FIG. 111-39 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR METHANE - FLOX PROPELLANT SYSTEM 

CH 4 - FLOX (82.6% F 2 ) 

P c = 300 PSIA (2.069 X 10 6 N/m 2 ) 

0/F = 6.50 
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FIG. 111-40 
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FIG. 111-41 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR METHANE - FLOX PROPELLANT SYSTEM 


CH - FLOX (82.6% F„ ) 

4 2 

P c = 500 PSIA (3.448 X 10 6 N/m 2 ) 
0/F = 5.50 



VACUUM SPECIFIC IMPULSE, T SP - (N) (SEC) / (k g) 

if VAC 
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FIG. 111-42 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR METHANE - FLOX PROPELLANT SYSTEM 

CH - FLOX (816% F, ) 

4 i 

P c « 500 PSIA (3.448 X 10 6 N/m 2 } 

O/F = 6.50 


(a/a^in) exit 
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FIG. ill -43 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
* FOR METHANE - FLOX PROPELLANT SYSTEM 


CH - FLOX (82.6% F, ) 

4 l 

750 PSIA (5.171 X 10 6 N/m 2 ) 
0/F = 3.50 
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FIG. 111-44 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR METHANE - FLOX PROPELLANT SYSTEM 


CH - FLOX (82.6% F„ ) 

4 2 

P c = 750 PSIA (5.171 X 10 6 N/m 2 ) 
0/F = 5.50 
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FIG. 111-45 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR METHANE - FLOX PROPELLANT SYSTEM 

CH 4 - FLOX (92.6% ) 

P c = 750 PSIA (5.171 X 10 6 N/m 2 ) 

0/F = 6.50 


o> 

-X 


U 

Ul 


Z 

I 


u 

< 

> 

0- 


m 

to 


3 

Q. 

U 

LL 

U 

Ul 

o. 


U 

> 




H9 10666- 16 


FIG. 111-46 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR METHANE - FLOX PROPELLANT SYSTEM 

CH, - FLOX (82.6% F, ) 

4 2 

P c =1000 PSIA (6.895 X 10 6 N/m 2 ) 

0/F = 4.00 
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FIG. 111-47 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR METHANE - FLOX PROPELLANT SYSTEM • 

CH 4 - FLOX (82.6% F 2 ) 

P c =1000 PSIA, (6.895 X 10 6 N/m 2 ) 

0/F = 6.00 
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FIG. 111-48 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR METHANE - FLOX PROPELLANT SYSTEM! 

CH 4 - FLOX (82.6% ? 2 ) 

P c =1000 PSIA (6.895 X 10 6 N/m 2 ) 

0/F = 6.50 
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TABLE IV - 1 

SUMMARY OF ELEMENTARY REACTIONS AND REACTION RATE CONSTANTS 
EMPLOYED IN AEROZINE 50-N 2 0 4 RECOMBINATION MECHANISM 


REACTION 

FORWARD RATE 

H + H + M H 2 +M 

Kf = 4.62 X lO 14 T _I 

1 h + h + h h 2 + H 

Kf =2 (4.62 X 10 U T" 1 ) 

2 H + H + H 2 h 2 +H 2 

Kf = 2.5 (4.62 X 10 14 T -1 ) 

3 H + H4-H 2 0 - 7 -— h 2 +h 3 o 

Kf = 2.5 (4.62 X 10 14 T -1) 

4 H + H + C0 2 — H 2 +C0 2 

Kf =2(4.62 X 10 14 T- 1 ) 

5 H + H + CO H 2 + CO 

Kf =2 (4.62 X 10 14 T- 1 ) 

6 H + H + N 2 H 2 + N 2 

Kf =2 (4.62 X 10 14 T-l) 

H + OH + M H 2 0 + M 

Kf = 7.85 X 10 ’S T -1 

7 h + OH+H — — h 2 o + h 

K f =3(7.85 X 10 ,s T" 1 ) 

8 H-fOH + H 2 H 2 0 + H 2 0 

Kf =3 (7.85 X 10 15 T _1 ) 

9 H +OH + H a O ^ 7 — HjO + HjO 

Kf = 20 (7.85 X 10 15 T _1 ) 

10 h + oh + co 2 ^^t h 2 o+co 

K f = 10 (7.85 X 1015 T" 1 ) 

11 H+OH + CO — H 2 0 + C0 

Kf =3 (7.85 X 10 1S T _1 ) 

12 H + OH + N 2 h 2 o + n 2 

Kf = 3 (7.85 X 10 15 T _1 ) 


NOTE: ALL RATES ARE EXPRESSED IN TERMS OF LB-MOLES, FT 3 , SEC; AND T IS IN °R 
TO CONVERT TO kg-MOLES, m 3 , SEC, AND T IN °K; MULTIPLY ABOVE RATE BY 2. 165 X 10~ 3 
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TABLE TV-2 

EQUILIBRIUM AND FROZEN FLOW PROPERTIES AT 
NOZZLE THROAT 

AEROZINE 50 (i) - N 2 0^(i) 



P^-psia 

OXIDIZER - FUEL WEIGHT RATIO 

1.25 

1.50 

1.75 

2.00 

2.25 

2.50 

Equilibrium 

Mass 

Flow - 

W/A 

Lom 

Ft 2 . Sec 

100 

81.8 

81.0 

81.5 

82.6 

84,1 

85.7 

200 

163.4 

161.5 

161.9 

164.1 

167.0 

170.2 

300 

245.0 

241.8 

242.1 

245.2 

249.5 

25^.3 

!■ 

4o8.1 

402.2 

402.1 

4o6.6 

■SH 

421.8 

750 

612.0 

602.6 

601.5 

607.8 

618.4 

630.5 

1000 

8.15.8 

802.8 

800.6 


■m 

828.6 








Frozen 

Mass 

- 

Lom 

Ft 2 . Sec 

100 

83.1 

83.0 

HI 

85.2 

HUH 

88.4 

200 

165.5 

164.9 

166.3 

168.9 

172.0 

175.2 

300 

mm 

246.4 


WBBKM 

B8S89I 


O 

O 

Lf\ 

411.9 

409.0 

411.5 

417.4 

425.0 

msm 

750 

616.8 

611.6 

6l4.6 

623.1 

634.4 

646.7 

1000 

821.6 

813.8 

817.1 

828.1 

843.0 

859.5 








Frozen 

Specific 

Heat 

Ratio 

>^R 

100 

BPS5S1 

1.24 

■m 

■PSi 

mmm 

mm 

200 

1.25 

1.24 

1.23 

1.23 

1.23 


300 

1.25 

1.24 

1.23 . 

■RSI 

BSH 

1BBS1 

500 

1.25 

1.24 

1.23 

1.23 

1.22 
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NOTE: P c (psia) x 6.895 x Io3 = P c (N/m 2 ) 

Mass Flow ( —2^ — ) x 4.883 = Mass Flow ( — j£S ) 

FtP.sec m2 sec ' 

















































































































H9 10666- 16 F | G> 

VARIATION OF EQUILIBRIUM VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

AEROZINE 50 (£) - N 2 0 4 (J>) 

P c = 100 PSIA <6.895X105N/m 2) 



EXIT AREA RATIO, (A/A Mm ) EX)T 
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FIG. IV-2 


VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

AEROZINE 50 (l) - N 2 0 4 (i) 

P c =100 PSIA(6,895X10 5 N/m2) 
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VARIATION OF EQUILIBRIUM VACUUM SPECIFIC IMPULSE WITH 

AEROZINE 50 (£) - N 2 0 4 (£) . 

P c = 200 PSIA (1.379X10^H/m 2) 
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H9 10666- 16 FIG, !V-5 

VARIATION OF EQUILIBRIUM VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

AEROZINE 50 (i) - N 2 0 4 (£) 

P c =300 PSIA (2.069X106N/m 2) 
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6 


VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

AEROZINE 50 {£) - N 2 0 4 (i) 

P c = 300 PSIA (2.069X106N/m2) 


























H9 10666- 16 F\&» |V-7 

VARIATION OF EQUILIBRIUM VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

AEROZINE 50 (£) - N 2 0 4 (£) 

P c = 500 PS!A (3.448X10 6 N/m 2) 
















H9 10666- 16 pig. IV-8 

VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

AEROZINE 50 (£) - N 2 0 4 (j?) 

P c = 500 PSIA (3.448X106H/m2) 
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H9 10666- 16 FIG. IV- 11 

VARIATION OF EQUILIBRIUM VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

AEROZINE 50 (£) - N 2 0 4 (I) 

P C =1000PSIA (6.895X106N/m 2) 
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H9 10666- 16 FIG, IV-12 

VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

AEROZINE 50 (i) - N 2 0 4 (i) 

P c =1000 PSIA (6.895X10 6 N/m?) 



EXIT AREA RATIO, (A/A M1N ) EX(T 

I 


















H ATOM GRADIENT, 0[H] /3t) - kg-MOLE m~ 3 SEC - 
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FIG. IV— 13 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

AEROZINE 50 (£) - N 2 0 4 (i) 

P c = 100 PSIA (6.895 X 10 5 N/m 2 ) 

O/F = 1.50 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE IV-1 



H ATOM GRADIENT, 0[Hj /8t) - kg-MOLE m" 3 SEC” 
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FIG. IV— 14 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

AEROZINE 50 (i) - n 2 o 4 (l) 

P c = 100 PSIA (6.895 X 10 5 N/m 2 ) 

O/F = 2.00 


NOTEs REACTION RATE CONSTANTS LISTED IN TABLE IV-1 





H ATOM GRADIENT, 0[H] /dt) - kg-MOLE m“ 3 SEC - 
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FIG. IV-15 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

AEROZINE 50 (i ) - N 2 0 4 (£) 

P c = 100 PSIA (6.895 X 10 5 N/m 2 ) 

O/F = 2.50 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE IV— 1 





H9 10666- 16 


FIG. IV — 16 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

AEROZINE 50 (1) - N 2 0 4 (£) 

P C = 200 PSIA 0.379 X 10 6 N/m 2 ) 

0/F = 1.50 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE IV— 1 
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FIG. IV— 17 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

AEROZINE 50 (£) - N 2 0 4 (i) 

P c = 200 PSIA (1.379 X 10 15 N/m 2 ) 

0/F = 2.00 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE IV-1 





H ATOM GRADIENT, 0[H] /at) - kg-MOLE m" 3 SEC" 


H9 10666- I 6 


FIG. IV-18 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

AEROZIN E 50 (£) - N 2 0 4 (I) 

P c = 200 PSIA (1.379 X 10 6 N/m 2 ) 

0/F = 2.50 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE IV-1 




H ATOM GRADIENT, (9[h] /9t) - kg-MOLE m~ 3 SEC" 
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FIG. IV— 19 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

AEROZINE SO (i) - N 2 0 4 (l) 

P c = 300 PSIA (2.069 X 10 6 N/m 2 ) 

O/F = 1.50 


NOTE; REACTION RATE CONSTANTS LISTED IN TABLE IV-1 
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FIG. IY-20 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

AEROZINE 50 (i) - N 2 0 4 (£) 

P c = 300 PSIA (2.069 X I0 6 N/m 2 ) 

0/F = 2.00 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE 1V-1 




H910666-16 FIG. IV— 21 

NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

AEROZINE SO(£) - N 2 0 4 (i) 

P c = 300 PSIA (2.069 X 10 6 N/m 2 ) 

0/F = 2.50 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE IV- 1 
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FIG. IV— 22 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

AEROZINE 50 (i) - N 2 0 4 (I) 

P c = 500 P SI A (3.448 X 10 6 N/m 2 ) 

0/F = 1.50 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE IV-1 































H ATOM GRADIENT, 0[H] /3t) - kg-MOLE m~ 3 SECT 
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FIG. IV-23 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

AEROZINE 50 (jO -N 2 0 4 (£) 

P c = 500 P SI A (3.448 X 10 6 N/m 2 ) 

O/F = 2.00 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE IV— 1 
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FIG. IV -24 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

AEROZINE 50 (i) - N 2 0 4 (£) 

P c = 500 PSIA (3.448 X 10 6 N/m 2 ) 

0/F = 2.50 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE IV-1 




H9 10666- 16 


FIG. IV —25 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

AEROZINE 50 (i) - N 2 0 4 (f) 

P c = 750 PSIA (5.171 X 10 6 N/m 2 ) 

0/F = 1.50 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE IV- 1 
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FIG. IV-26 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

AEROZINE 50 (1) - N 2 0 4 {£) 

P c = 750PSIA (5.171 X 10 6 N/m 2 ) 

0/F = 2.00 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE IV— I 
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FIG. IV— 27 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

AEROZINE 50 (£) - N 2 0 4 (i) 

P c =750 PSIA (5.171 X 10 6 N/m 2 ) 

O/F =2.50 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE IV— 1 


U 

US 

VO 


ro 

£ 

ULi 


o 

I 


cn 

JK 


I 


CO 

\ 


3E 

CO 


z 

IU 

a 

< 

oc 

o 

o 




H ATOM GRADIENT, 0[H] /dt) - kg-MOLE nr 3 SEC" 


H9 10666- 16 


FIG. IV-28 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

AEROZINE S0(£) - N 2 0 4 (£) 

P c =1000 PSIA (6.895 X 10 6 N/m 2 ) 

0/F = 1.50 
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FIG. IV-29 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

AEROZINE 50(|) -N 2 0 4 (i) 

P C =1000PSIA (6.895 X 10 6 N/m 2 ) 

0/F = 2.00 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE IV— 1 
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FIG. IV-30 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

AEROZINE 50 (i) - N 2 0 4 (i) 

P c =1000 PSI A (6.895 X 10 6 N/m 2 ) 

0/F =2.50 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE IV- 1 
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FIG. IV— 31 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR AER0Z1NE -N 2 0 4 PROPELLANT SYSTEM 

AEROZINE 50 <i) - Nj 0 4 (i ) 

P c =100 PSIA (6.895X105N/m 2) 

0/F = 1.50 



FREEZING AREA RATIO, (A/A M1N ) 
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FIG. IV-32 


EFFECT OF FREEZING AREA RATIO ON NONEQUIUBRIUM PERFORMANCE 
FOR AEROZINE - N2 0 4 PROPELLANT SYSTEM 


AEROZINE 50 (I) - N 2 0 4 (i) 
P c =100 PSIA <6.895X105N/m2) 
0/F = 2.00 
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FIG. IV— 33 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR AEROZINE -N 2 0 4 PROPELLANT SYSTEM 


AEROZINE 50 (i) - N 2 0 4 (I) 
P c = 100 PSIA (6.895X105N/m 2 ) 
0/F - 2.50 
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FIG. IV-34 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR AEROZINE -N 2 0 4 PROPELLANT SYSTEM 


AEROZINE 50 (£) - N 2 0 4 (i ) 
P c = 200 PSIA (1.379X106N/m2) 
0/F = 1.50 
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FIG. IV-35 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR AEROZINE -N 2 0 4 PROPELLANT SYSTEM 


AEROZINE 50 (£) - N 2 0 4 (£) 
P c = 200 PSIA (1.379X10 6 N/m 2 ) 
0/F = 2.00 
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FIG, IV— 36 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR AEROZINE -N 2 0 4 PROPELLANT SYSTEM 

AEROZINE 50 (i) - Nj 0 4 (f) 

P c = 200 PSIA (1.379X10 6 N/m 2 ) 

0/F = 2.50 
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FIG. IV - 37 


EFFECT OF FREEZING AREA RATIO ON NONEQUILf BRIUM PERFORMANCE 
FOR AEROZINE -N 2 0 4 PROPELLANT SYSTEM 


AEROZINE 50(f) - N 2 0 4 (f) 
P c = 300 PSIA ( 2.069 X 10 6 N /m 2) 
0/F = 1.50 
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FIG. IV -38 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR AEROZINE -N2O4 PROPELLANT SYSTEM 


AEROZINE 50 (i) - N 2 0 4 (i) 
P c = 300 PSIA ( 2 . 069 X 10 <Wm 2 ) 
0 /F = 2.00 
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FIG. IV-39 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR AEROZINE -N 2 0 4 PROPELLANT SYSTEM 

AEROZINE 50 (i) - Nj 0 4 (I) 

P c = 300 PSIA ( 2.069 X 10 6 N/m 2 ) 

0/F = 2.50 
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FIG. IV-40 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR AEROZINE -N 2 0 4 PROPELLANT SYSTEM 

AEROZINE 50(i) - N 2 0 4 (£) 

P c = 500 PSIA (3.448X10 6 N/m 2 ) 

O/F = 1.50 


360 




i 


u 

< 


> 340 

Q. 





u 

lu 

<l 320 

i 

3 

U 

< 

> 


310 


300 



FREEZING AREA RATIO, (A/A MIN > 

r R 









VACUUM SPECIFIC IMPULSE, I, p - (N) (SEC)/ (Vg) 

bK VAC 


H9 10666- 16 


FIG. IV— 4 1 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR AEROZINE -N 2 0 4 PROPELLANT SYSTEM 

AEROZINE 50 <£) - N 2 0 4 <f ) 

P c = 500 PSIA (3.448X10 6 N/m2) 

0/F = 2.00 
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FIG. IV -42 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR AEROZINE -N 2 0 4 PROPELLANT SYSTEM 

AEROZINE 50 (i) - N 2 0 4 (ft) 

P C = 500 PSIA (3.448X10 6 N/m 2 ) 

0/F = 2.50 
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FIG. IV-43 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR AEROZINE -N 2 0 4 PROPELLANT SYSTEM 

AER0Z1NE 50 (£) - Nj 0 4 (i ) 

P c = 750 PSIA (5.171X10 6 N/m 2 ) 

0/F = 1.50 
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FIG. IV-44 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR AEROZINE -N 2 0 4 PROPELLANT SYSTEM 

AEROZINE 50 (t ) - N 2 0 4 (i) 

P c = 750 PSIA (5.171X10 6 N/m 2 ) 

O/F = 2.00 
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FIG. IV-45 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR AEROZINE -N 2 0 4 PROPELLANT SYSTEM 

AEROZINE 50 (I) - N 2 0 4 (£ ) 

P c = 750 PS!A(5.171X10 6 N/m 2 ) 

0/F = 2.50 
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FIG. IV -46 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR AEROZINE - N 2 0 4 PROPELLANT SYSTEM 


AEROZINE 50 (i ) - N 2 <3, (I) 
P c =1000 PSIA (6.895X10 2 3 * * 6 N/m2) 
0/F = :.50 
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FIG. IV-47 


\ 

EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR AEROZINE -N2O4 PROPELLANT SYSTEM 


AEROZINE 50 {£) - N 2 0 4 (£) 
P c =1000 PSIA (6.895X10 6 N/m2) 
0/F = 2.00 
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FIG. IV-48 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR AEROZINE -N 2 0 4 PROPELLANT SYSTEM 

AEROZINE 50 (2) - N 2 0 4 (£) 

P c =1000 PSIA (6.895X10 6 N/m 2) 

0/F - 2.50 
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TABLE 2-1 

SUMMARY OF ELEMENTARY REACTIONS AND REACTION RATE CONSTANTS 
EMPLOYED IN B^-OF 2 RECOMBINATION MECHANISM 


REACTION 

FORWARD RATE 


H + H + A 


H^f A 

k f = 4.62 X 10 14 T-1 

1 

H + H + H 


Hj + H 

k f = 2.0 (4.62 X 10 14 T- 1 ) 

2 

H + H + H 2 


h 2+ h 2 

kf = 2.5 (4.62 X I0 14 T- 1 ) 

3 

H + H + HF 

— 

h 2 + hf 

k f = 2.5 (4.62 X 10 14 T” 1 ) 

4 

H + H + BOF 


Hj + BOF 

k f = 2.0(4.62 X 10’ 4 T“>) 

5 

H + H + BF 

— 

H 2 + BF 

kf = 2.0 (4.62 X 10 14 T” 1 ) 

6 

H + H + H^O 


^ +HJO 

k, = 2.5 (4.62 X 10 14 T- 1 > 


H + F + A 


HF + A 

kf = 1.155 X 10 15 T" 1 

7 

H + F + H 

— 

HF + H 

k f = 1 (1.155 X 10 15 T-> ) 

8 

H + F + H 2 


HF .+ H 2 

k f = - 2.5 (1.155 X 10 15 T- 1 ) 

9 

H + F + HF 

— 

HF 4- HF 

kf a 2.5 (1.155 X 10' S T - 1 ) 

10 

H + F + BOF 


HF + BOF 

k,= 2.0 (1.155 X 10 15 T- 1 ) 

11 

H + F + BF 

— 

HF + BF 

kf = 2.0 (1.155 X 10 I s T** 1 ) 

12 

H + F + HjO 


HF + ^0 

kf = 2.5 (1.155 X 10 15 T“1 ) 


H + OH + A 


H 2 0 + A 

kf = 7.85 X 10 15 T“1 

13 

H +OH + H 

-*• 

HjO + H 

kf = 3 (7.85 X 10 15 T“ 1 ) 

14 

H + OH + H 2 

— 

h 2 o+ h 2 

kf = 3 (7.85 X 10 15 T _1 ) 

15 

H 4- OH +HF 


H 2 0 + HF 

kf = 3 (7.85 X 10 15 T“' ) 

16 

H + OH + BOF 

— 

h 2 o + bof 

kf = 10 (7.85 X 10 15 T" 1 ) 

17 

H + OH + BF 

— 

H 2 0 + BF 

kf = 3 (7.85 X 10 ,s T“' ) 

18 

h + oh + h 2 o 

— 

h 2 o + h 2 o 

kf = 20 (7.85 X 10 ,s T“* ) 


NOTE: ALU RATES EXP RESSED IN TERMS 0 F LB-MOLES, FT 3 , SEC; AND T IS IN ° R 


TO CONVERT TO Kg-MOLES, M 3 . SEC. AND T IN °K. MULTIPLY ABOVE RATE BY 2.165 X 10“ 3 
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TABLE V-2 


EQUILIBRIUM AND FROZEN PROPERTIES AT 
NOZZLE THROAT 
B 2 H 6 (i) - OF 2 ( ^ ) 



P c -psia 

OXIDIZER - FUEL WEIGHT RATIO 

2.5 

3.0 

3.5 

4.0 

4.5 

5.0 

5.5 

Equilibrium 
Mass 
Flow -• 
Lbm 

50 

33 « 5 ^ 

33.16 

32.92 

32.90 

33.07 

33.33 

33.61 

100 

66.46 


65.20 

65.16 

65.49 

65.98 

-- . 

250 

164.2 

162.2 

161.0 

160.9 

161.6 

162.7 

164.1 

500 

325 . T 

321.4 

319.1 

318.7 

■EaglW 

322.3 

324.8 

Ft^.Sec 









Frozen 
Mass 
Flow - 
Lbm 

50 

urea 

wmam 

Wifi 

34.47 

nisa 

35.03 

" 35 . 3 ^ 

100 

68.98 

68.27 

67.96 

68.11 

68.57 

BEgEB 

-- 

250 

169.9 

168.0 

167.2 

167.5 

168.6 

170.0 

171.6 

500 

336.1 

332.2 

330.6 

■asw 

333.1 

335.7 

338.7 

Ft 2 .Sec 









Frozen 
Specific i 
Heat 
Ratio 

50 

1.278 

1.286 

1.292 


1.300 

1.303 

1.303 

100 

1.274 

1.281 

1.287 

1.291 


1.298 

— 

250 

1.264 

1.271 

1.276 

1.281 

1.284 

1.288 

1.291 

500 

1.263 

1.269 

1.274 

1.278 

1.282 

1.282 

1.288 


NOTE: P c (psia) x 6.895 x 10 ^ = P c (N/m^) 

Mass Flow ( -- w— ) x 4.883 = Mass Flow 

\Ft 2 .sec ' 



kg 

sec 


) 
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FIG. V-1 


VARIATION OF EQUILIBRIUM VACUUM SPECIFIC IMPULSE WITH AREA RATIO 


B 2 h 6 (j?) - 0F 2 (i) 

P c =50 PSIA (3.448 X 10 5 N/m 2 ) 
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FIG. V-2 


VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 


B 2 H 6 (i) - 0F 2 (£) 

P c = 50 PSIA (3.448 X 10 5 N/m 2 ) 




VACUUM SPECIFIC IMPULSE, l SP - (N) (SEC)/(kg) 
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FIG. V-3 
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FIG. V-4 


VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 


b 2 h 6 a) - of 2 (f> 

P c =100 PSIA (6.895 X 10 5 N/m 2 ) 



EXIT AREA RATIO, (A/A M1N ) EX1T 


VACUUM SPECIFIC IMPULSE, I S p ( vac) - (N)(SEC)/(kg) 


H910666-16 


FIG. V-5 


VARIATION OF EQUILIBRIUM VACUUM SPECIFIC IMPULSE WITH AREA RATIO 


B 2 H 6 (£) - OF 2 (£) 

P c = 250 PSIA (1.724 X 10 6 N/m 2 ) 









VACUUM SPECIFIC IMPULSE, Jc P/ „ A _, - (N) (SEC)/(kg) 

J ( V AC ) 
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FIG. V-6 


VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 

b 2 h 6 m - of 2 (f> 

Pc = 250 PSIA (1.724 X IQ 6 N/tn 2 ) 



20 


80 
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FIG. V-7 


VARIATION OF EQUILIBRIUM VACUUM SPECIFIC IMPULSE WITH AREA RATIO 


B 2 H 6 U) - 0F 2 (I) 

P c =500 PSIA (3.448 X 10 6 N/.m 2 ) 



i 




VACUUM SPECIFIC IMPULSE, l S p (VAC) ~ (N)(SEC)/(kg) 
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FIG. V-8 


VARIATION OF FROZEN VACUUM SPECIFIC IMPULSE WITH AREA RATIO 


B 2 H 6 a) - OF 2 (£) 

■P c =500 PSIA (3.488 X 10 6 N/nr 2 ) 



EXIT AREA RATIO, (A/A M1N ) EX1T 




H ATOM GRADIENT ( d [H]/<9t) p - kg - MOLE m SEC 
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FIG. V-9 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 

B 2 H 6 (£) -OF 2 (I) 

P c = 50 PSIA (3.448 X 10 5 N/m 2 ) 

O/F = 3.0 

NOTF.: REACTION RATE CONSTANTS LISTED IN TABLE V-l 




H ATOM GRADIENT (d [H]/dt) p- kg - MOLE m SEC ~ 


H9 10666- 16 


FIG. V-10 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


b 2 h 6 U) - of 2 u) 

P c = 50 PSIA (3.448 X 10 5 N/m 2 ) 


0/F = 3.5 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE V-l 




H ATOM GRADIENT (d [H]/<9t) p - kg - MOLE m SEC 
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FIG. V-11 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


B 2 H 6 (£) -OF 2 (£) 

P c = 50 PSIA (3.448 X 10 5 N/m 2 ) 


0/F = 4.0 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE V-l 




H ATOM GRADIENT ( d [H]/d*) p - kg - MOLE m SEC “ 
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FIG. V-12 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


B 2 H 6 (?) - OF 2 (1) 

P c = 50 PSIA (3.448 X 10 5 N/m 2 ) 
0/F = 4,5 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE V_1 




H ATOM GRADIENT (<9 [H]/<9t) p- kg - MOLE m~ J SEC* 
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FIG. V-13 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


B 2 h 6 (£) - 0F 2 (i) 

P c - 100 PS1A (3.448 X 10 5 N/m 2 ) 
0/F = 3.0 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE V-l 






H ATOM GRADIENT (d [H]/c?t) p - kg - MOLE m SEC 
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FIG. V-14 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


B 2 H 6 ( 1 ) - O F 2 <■*> 

P c = 100 PS1A (3.448 X 10 5 N/m 2 ) 
0/F = 3.5 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE V-l 




H ATOM GRADIENT ( d [H]/3t) p- kg - MOLE m SEC' 
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FIG. V-15 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


b 2 h 6 (j?) - of 2 a) 

P c = 100 PSIA (3.448 X 10 5 N/m 2 ) 
O/F = 4.0 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE V-l 





H ATOM GRADIENT ( d [H]/<?t) p - kg - MOLE m SEC ' 
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FIG. V-16 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


B 2 H 6 (i) -OF 2 «) 

P c -100 PS1A (3.448 X 10 5 N/m 2 ) 
0/F =4.5 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE V-l 





H ATOM GRADIENT (<5 [H]/<Jt) p - kg- MOLE m SEC * 
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FIG. V-17 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


B 2 H 6 (£) - OF 2 (i) 

P c = 250 PSIA (3.448 X 10 5 N/m 2 ) 


0/F = 3.0 




H ATOM GRADIENT (d [HVdl) p - kg - MOLE m SEC ~ 
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FIG. V-18 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


b 2 h 6 (i) - 0F 2 (j?) 

P c - 250 PSIA (3.448 X 10 5 N/m 2 ) 


O/F = 3.5 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE V- 1 




H ATOM GRADIENT (d [H]/<9t) p - kg - MOLE SEC' 
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FIG. V-19 


10 3 


- 10 2 


10 * 


10 ° 


10 “ 1 


10 “ 2i 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


B 2 H 6 (1) - OF 2 (t) 

P c = 250 PSIA (3.448 X 10 5 N/m 2 ) 
0/F = 4.0 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE V-1 



H ATOM GRADIENT ( d [H]/<9t) p- kg - MOLE m SEC 
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FIG. V-20 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


B 2 H 6 (j?) - OF 2 (I!) 

P c = 250 PSIA (3.448 X 10 5 N/m 2 ) 
O/F =4.5 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE V-l 




H ATOM GRADIENT ( d [H]/dt) p- kg - MOLE m SEC ~ 
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FIG. V-21 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


B 2 H 6 (D - OF 2 (O 
P c = 500 PSIA (3.448 X 10 5 N/m 2 ) 
O/F = 3.0 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE V-l 





H ATOM GRADIENT (d [H]/£t) p- kg - MOLE m " SEC ~ 
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FIG. V-22 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


B 2 H 6 U> -OF 2 (i) 

P c =. 500 PSIA (3.448 X 10 5 N/m 2 ) 
0/F = 3.5 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE V-l 




H ATOM GRADIENT (d [HVdi) p- kg - MOLE m SEC “ 
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FIG. V-23 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


B 2 H 6 (jH - OF 2 ti) 

P c = 500 PSIA (3.448 X 10 5 N/m 2 ) 
0/F = 4.0 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE V-l 





H ATOM GRADIENT (<9 [H]/<9t) p- kg - MOLE SEC ' 
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FIG. V-24 


NORMALIZED GRAPHICAL SOLUTION FOR FREEZING AREA RATIO 
USING MODIFIED BRAY ANALYSIS 


B 2 H 6 (£) - OF 2 (f) 

P c = 500 PSIA (3.448 X 10 5 N/m 2 ) 
0/F = 4.5 


NOTE: REACTION RATE CONSTANTS LISTED IN TABLE V-l 
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FIG. V-25 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR A DIBORANE - 0F 2 PROPELLANT SYSTEM 

B 2 H 6 ({) - of 2 (f) 

P c = 50 PSIA (3.448 X 10 5 N/m 2 ) 

0/F =3.0 
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FIG. V-26 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR A DIBORANE - OF 2 PROPELLANT SYSTEM 

B 2 h 6 (£) - of 2 (f) 

P c = 50 PSIA (3.448 X 10 5 N/m 2 ) 

0/F = 3.5 



FREEZING AREA RATIO (A/A M1N ) FR 


100 

80 

60 


40 


20 





H910666-16 


FIG. V-27 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR A DIBORANE - OF 2 PROPELLANT SYSTEM 

B 2 H 6 (£) - OF 2 (£) 

P c = 50 PSIA (3.448 X IQ 5 N/m 2 ) 

0/F = 4.0 
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FIG. V-28 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR A DIBORANE - 0F 2 PROPELLANT SYSTEM 

B 2 H 6 (f) - 0F 2 ( 2 ) 

P c = 50 PSIA (3.448 X 10 5 N/m 2 ) 

0/F =4.5 
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FIG. V-29 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR A DIBORANE - 0F 2 PROPELLANT SYSTEM 

B 2 H 6 <£) - 0F 2 (I) 

P c = 100 PSIA ( 6.985 X 10 5 N/m 2 ) 

0/F =3.0 
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FIG. V-30 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR A DIBORANE - 0F 2 PROPELLANT SYSTEM 

b 2 h 6 <f) - 0F 2 (() 

P c = 100 PS1A ( 6.985 X 10 5 N/m 2 ) 


0/F = 3.5 
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FIG. V-31 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR A DIBORANE - 0F 2 PROPELLANT SYSTEM 

B 2 H 6 (£) - OF 2 (f) 

P c = 100 PSIA ( 6.985 X 10 s N/m 2 ) 

0/F = 4.0 
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VACUUM SPECIFIC IMPULSE, l S p (VA c) ” (N) < S EC)/<k 9 ) 
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FIG. V-32 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR A DIBORANE - 0F 2 PROPELLANT SYSTEM 

B 2 H 6 (l) - of 2 ({') 

P c = 100 PSIA ( 6.985 X 10 s N/m 2 ) 

O/F =4.5 





FREEZING AREA RATIO (A/A M1N ) FR 





VACUUM SPECIFIC IMPULSE, »SP( VAC ) " M (SEC)/(kg) 
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FIG. V-33 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR A DIBORANE - 0F 2 PROPELLANT SYSTEM 

b 2 h 6 ({) - of 2 (i) 

P c = 250 PSIA ( 1.742 X 10 6 N/m 2 ) 

0/F =3.0 



FREEZING AREA RATIO (A/A M|N ) FR 






VACUUM SPECIFIC IMPULSE, ls P(VAC) - (N) (SEC)/(kg) 
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FIG. V-34 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR A DIBORANE - 0F 2 PROPELLANT SYSTEM 

B 2 H 6 (£) - 0F 2 ({) 

P c = 250 PSIA ( 1.742 X 10 6 N/m 2 ) 

O/F = 3.5 





FREEZING AREA RATIO (A/A M|N ) FR 





VACUUM SPECIFIC IMPULSE, Isp ( vac) ~ ( N > (SEC)/(kg) 
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FIG. V-35 


44001 


4300 


4200 


4100 


4000 


3900 


3800 


3700 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR A DIBORANE - 0F 2 PROPELLANT SYSTEM 

B 2 H 6 (£) - 0F 2 (£) 

P c = 250 PSIA ( 1.742 X 10 6 N/m 2 ) 

O/F =4.0 
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FREEZING AREA RATIO (A/A M1N ) FR 
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FIG. V-36 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR A DIBORANE - OF 2 PROPELLANT SYSTEM 

B 2 H 6 (£) - OF 2 (£) 

P c = 250 PSIA ( T.742 X 10 6 N/m 2 ) 

0/F = 4.5 
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FREEZING AREA RATIO (A/A MIN ) FR 



VACUUM SPECIFIC IMPULSE, I SP (VAC) - (N) (SEC)/(kg) 
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FIG. V-37 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR A DIBORANE - OF 2 PROPELLANT SYSTEM 

B 2 H 6 <£) - OF 2 (I) 

P c = 500 PSIA (3.4475 X 10 6 N/m 2 ) 

O/F =3.0 





FREEZING AREA RATIO (A/A MIN ) FR 
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FIG. V-38 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR A DIBORANE - 0F 2 PROPELLANT SYSTEM 

B 2 H 6 (I) - 0F 2 «) 

P c a 500 PSIA (3.4475 X 10 6 N/m 2 ) 

0/F =3.5 
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FIG. V-39 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR A DIBORANE - 0F 2 PROPELLANT SYSTEM 

B 2 H 6 <£) - 0F 2 (i) 

P c = 500 PSIA (3.4475 X 10 6 N/m 2 ) 
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FIG. V-40 


EFFECT OF FREEZING AREA RATIO ON NONEQUILIBRIUM PERFORMANCE 
FOR A DIBORANE - 0F 2 PROPELLANT SYSTEM 

B 2 H 6 (£) - 0 F 2 (I) 

P c = 500 PS I A (3.4475 X 10 6 N/m 2 ) 

0/F =4.5 
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FIG. AL-1 
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FIG. AL-2 
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VARIATION OF DIVERGENCE LOSS AS PERCENT OF FROZEN 
VACUUM SPECIFIC IMPULSE FOR 15 DEG CONICAL NOZZLE 
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FIG. AL-4 



IMSDHSd - $$01 33N30H3AIG 


EXIT AREA RATIO, (A/A MIN 


VARIATION OF DIVERGENCE LOSS AS PERCENT OF FROZEN 
VACUUM SPECIFIC IMPULSE FOR 70% BELL NOZZLE 
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FIG. AL-5 










VARIATION OF VISCOUS DRAG AS PERCENT OF FROZEN VACUUM SPECIFIC IMPULSE 

FOR 15° CONICAL NOZZLE 
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FIG. AL-6 
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VARIATION OF VISCOUS DRAG AS PERCENT OF FROZEN VACUUM SPECIFIC IMPULSE 
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CONICAL NOZZLE 


H910666-16 


FIG. AL-8 
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VARIATION OF VISCOUS DRAG AS PERCENT OF FROZEN VACUUM SPECIFIC IMPULSE 

FOR 70% BELL NOZZLE 
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FIG. AL-9 
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VARIATION OF VISCOUS DRAG AS PERCENT OF FROZEN VACUUM SPECIFIC IMPULSE 

FOR 70% BELL NOZZLE 


H9 10666- 16 


FIG. AL-10 
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VARIATION OF VISCOUS DRAG AS PERCENT OF FROZEN VACUUM SPECIFIC IMPULSE 
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FIG. AL-12 
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FIG. AL-13 





